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ABSTRACT 

Two  design  methods  for  multi-parameter  dynamic  systems 
are  proposed.   They  are  intended  to  eliminate  the  limita- 
tions and  disadvantages  of  the  existing  design  methods.  The 
powerful  mathematical  tools  of  optimal  control  theory  are 
applied  to  the  practical  design  problems  of  classical  con- 
•  trol. 

The  first  method  is  intended   for  linear  systems  only; 
the  design  problem  is  solved  in  the  s-domain,  by  finding 
"the  best  root  locations"  of  the  system's  characteristic  e- 
quation.   In  the  second  method,  the  design  problem  is  solved 
by  finding  "the  best  response"  of  the  system  in  the  time 
domain .   The  second  method  is  applicable  to  a  wide  range  of 
dynamic  systems;   it  can  be  used  to  synthesize  linear,  non- 
linear and  sampled-data  systems,  and  systems  with  time  de- 
lay.  This  method  is  also  extended  to  a  numerical  stability 
analysis  procedure. 

Fourteen  examples  are  presented  to  illustrate  the 
applications  of  the  methods. 
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I.   INTRODUCTION 

Classical  control  theory  and  modern  control  theory 
have  been  extended  in  diverging  routes  since  the  inception 
of  modern  control  theory.   Modern  techniques  have  been  suc- 
cessful for  larger  scale  problems  and  for  problems  which 
include  some  new  performance  criteria,  such  as  expenditure 
of  control  energy,  etc.,  but  a  general  purpose  or  a  special 
purpose  computer  in  the  control  loop  is  usually  necessary 
to  implement  the  solution. 

Most  design  techniques  of  classical  control  theory 
were  developed  before  the  extensive  application  of  digital 
computers  to  system  design  problems.   Using  a  digital  com- 
puter in  the  design  phase  of  a  control  system  is  usually 
feasible,  especially  if  the  design  can  be  implemented  by 
using  passive  circuit  elements  rather  than  complicated  and 
expensive  memory  units. 

During  the  last  decade  mathematical  methods  and  numer- 
ical teqhniques  have  been  developed  and  applied  to  the  prob- 
lems of  "Optimal  Control  Theory".   The  same  methods  have 
found  application  in  other  areas  of  science,  but  they  have 
not  been  applied  efficiently  to  the  design  procedures  of 
the  dynamic  systems  which  have  been  defined  in  the  sense  of 
classical  control  theory.   The  classical  techniques,  al- 
though having  many  drawbacks  and  limitations  (see  Section 
II. A),  were  well  developed  by  the  time  effective  optimiza- 
tion methods  arrived  on  the  scene.   The  classical  methods 


have  been  successful  for  a  restricted  class  of  design  prob- 
lems;  for  this  reason  probably,  modern  control  theory  grew 
independently,  rather  than  improving  and  replacing  the 
existing  methods. 

In  this  thesis  it  is  intended  to  narrow  the  gap  be- 
tween the  classical  and  modern  control  theories.   Powerful 
mathematical  tools  of  optimization  theory  are  applied  to 
the  practical  design  problems  of  classical  control.   Two 
general  design  methods  for  multi-parameter  dynamic  systems 
are  proposed: 

1.  Optimization  for  the  Best  Root  Location  in  the 

i 

s-Domain,  and 

2.  Optimization  for  the  Best  Response  in  the  Time 
Domain. 

Chapter  II  of  the  thesis  is  devoted  to  the  first  meth- 
od which  is  intended  for  linear  systems  only  —  the  design 
problem  is  solved  in  the  s-domain  by  finding  "the  best  root 
locations"  of  the  system's  characteristic  equation.   Four 
linear  system  design  examples  are  included,,   The  rest  of 
the  thesis  is  devoted  to  the  second  method  and  its  various 
applications,   In  Chapter  III,  a  general  description  and 
the  philosophy  of  the  method  "Optimization  for  the  Best  Re- 
sponse in  the  Time-Domain"  is  given,  and  it  is  applied  to 
the  same  linear  system  design  problems  discussed  in  Chapter 
II,  to  provide  a  basis  for  comparison  of  the  results  ob- 
tained by  both  methods. 
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The  time-domain  method  is  applicable  to  a  large  varie- 
ty of  system  design  problems.   To  demonstrate  this  point. 
Chapter  IV  is  devoted  to  various  applications  of  the  method. 
Three  examples  for  nonlinear  and  sampled-data  systems  and 
a  system  with  transport  lag  are  presented  in  Chapter  IV. 

In  Chapter  V,  a  numerical  method  to  investigate  the 
stability  limit  of  the  dynamic  systems.,  is  introduced  as  an 
extension  of  the  time-domain  method. 

In  the  year  1964 ,  the  gap  between  control  theory  and 
the  practical  control  problems  became  substantial  and  ex- 
tensive research  effort  has  been  directed  to  bridge  this 
gap  between  the  theory  and  paractice  (see,  for  example,  Ref. 
2) .   The  methods  presented  in  this  thesis  can  be  thought  of 
as  another  approach  to  suboptimal  system  design  problems. 
In  this  respect,  these  methods  can  be  considered  as  applica- 
tions and  extensions  of  Optimal  Control  Theory;  on  the  other 
hand,  these  methods  provide  better  solutions  to  the  practi- 
cal dynamic  system  synthesis  problems  of  Classical  Control 
and  they  can  be  considered  as  the  extensions  of  Classical 
Control.   Since  the  extensions  from  Optimal  Control  Theory 
and  Classical  Control  are  toward  each  other,  these  methods 
hopefully  provide  one  of  the  necessary  bridges  to  reduce 
the  gap. 


II.   LINEAR  SYSTEM  DESIGN  BY  OPTIMIZATION  FOR 
THE  BEST  ROOT  LOCATIONS  IN  THE  S-DOMAIN 

A.    GENERAL 

Present  design  methods  of  classical  control  theory  are 
limited  and  have  some  serious  disadvantages.   Root-locus 
and  frequency-domain  techniques,  for  example,  can  handle 
systems  which  have  only  one  parameter.   Algebraic  methods 
(in  the  coefficient  and  parameter  planes)  are  two-variable 
procedures.   The  major  design  tools  of  these  methods  are  a 
set(s)  of  graphs,  and  obtaining  these  graphs  for  higher- 
order  systems  is  quite  a  laborious  task.   The  trial-and- 
error  nature  of  the  classical  design  procedures  is,  however, 
the  most  significant  limitation  of  these  methods. 

Since  the  dynamic  behavior  of  second-order  systems  has 
been  thoroughly  studied  and  is  well  known,  most  of  the  pre- 
sent methods  for  high-order  system  design  depend  on  the 
dominant-roots  concepts.   Finding  a  pair  of  complex  dominant 
roots  for  a  high-order  system  is  equivalent  to  approximat- 
ing the  system  by  a  second-order  system. 

A  pair  of  complex  roots  of  high-order  system  can  be 
made  dominant  by  satisfying  at  least  one  of  the  two  condi- 
tions given  below: 

1.   Transients  in  the  time  response  due  to  undesired 
roots  (all  the  roots  except  the  selected  "dominant"  pair) 
last  for  a  much  shorter  time  than  the  transients  of  the  se- 
lected pair.   This  is  achieved  if  the  time  constants  of  the 
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undesired  roots  are  much  smaller  than  those  of  the  selected 
pair.   In  the  s-domain,  the  magnitudes  of  the  real  parts  of 
the  undesired  roots  must  be  much  larger  than  the  magnitu- 
des of  the  real  parts  of  the  selected  pair  to  satisfy  this 
condition.   The  remote  undesired  roots  may  cause  undesir- 
able effects  in  the  time  response  for  a  short  period  of 
time  (near  zero  time) ,  but  thereafter  the  time  response  is 
governed  by  the  selected  roots. 

2.  The  coefficients  in  the  time  response  associated 
with  the  desired  roots  are  much  larger  than  those  associ- 
ated with  the  undesired  roots.  In  other  words,  the  contri- 
bution of  the  selected  pair  to  the  transient  response  is 
much  larger  than  the  contribution  of  the  other  roots.  In 
the  s-domain,  this  condition  is  satisfied  if  the  residues 
of  the  desired  roots  are  much  larger  than  the  residues  of 
the  other  roots. 

If  both  of  these  conditions  are  satisfied,  the  time 
response  of  the  selected  pair  dominates  the  overall  time 
response,  and  the  selected  pair  of  roots  is  called 
"dominant".   If  only  one  of  the  conditions  is  satisfied  an 
acceptable  transient  response  may  still  be  obtained.   The 
first  condition  alone,  for  example,  has  been  used  in  the 
s-domain  design  procedures  to  avoid  the  lengthy  and  time- 
consuming  residue  calculations  at  each  step  of  the  trial- 
and-error  procedures. 

When  designing  a  high-order  feedback  system  by  root- 
locus  methods,  or  by  algebraic  methods,  the  designer 
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usually  selects  a  pair  of  complex  root  locations  in  the  s- 
domain,  and  then  tries  to  locate  two  of  the  roots  of  the 
system's  characteristic  equation  at  these  selected  loca- 
tions and  the  other  (undesired)  roots  at  remote  points  of 
the  left-hand  side  of  the  s-plane.   If  this  can  be  achieved, 
the  transient  response  of  the  system  closely  follows  the 
selected  (desired)  second-order  system's  response.   There 
are  usually  other  specifications  to  be  met,  such  as  the 
steady-state  error,  maximum  allowable  gain,  and  limitations 
on  the  system  components,  which  are  treated  as  additional 
constraints. 

The  number  of  free  independent  parameters  (N  )  cannot 
exceed  the  order  of  the  characteristic  equation  (n) ,  and  if 
the  number  of  free  independent  parameters  is  equal  to  the 
order  of  the  system  (N  =  n) ,  the  roots  of  the  characteris- 
tic equation  can  be  located  at  any  point  in  the  s-plane. 
The  design  is  completed  by  simple  algebraic  manipulations 
(see, for  example,  Ref.  12).   If  the  number  of  free  indepen- 
dent parameters  is  less  than  the  order  of  the  characteris- 
tic equation  of  the  system  (N   <  n) ,  however,  the  designer 
loses  his  control  over  some  of  the  roots.   Since  one  free 
parameter  exactly  locates  one  root,  when  all  parameters  are 
used  to  locate  an  equal  number  of  roots,  the  excess  (n  -  N  ) 
root(s)  may  be  anywhere  in  the  s-plane.   Some  of  these 
roots  may  move  very  close  to  the  origin  and  become  dominant, 
or  some  may  even  go  into  the  right  half  of  the  s-plane  and 
make  the  system  unstable. 
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Two  parameters  can  be  used  to  locate  a  pair  of  complex 
roots  exactly  at  the  specified  locations,  but  these  roots 
can  be  constrained  to  move  on  a  line  by  using  just  one  pa- 
rameter.  For  example,  their  real  parts  (cr)  or  imaginary 
parts  (a))  can  be  fixed,  or  they  can  be  forced  to  move  on  a 
constant  damping  ratio  (C)-line  or  on  a  constant  natural 
frequency  (oj  ) -circle. 

When  designing  with  a  fewer  number  of  parameters  than 
the  degree  of  the  characteristic  equation,  the  procedure  is 
a  trial-and-error  method  and  success  may  depend  primarily 
on  the  past  experience  of  the  designer.   There  are  several 
methods  for  locating  two  of  the  roots  of  the  characteristic 
equation  at  the  desired  locations,  but  if  the  number  of 
free  parameters  is  less  than  the  order  of  the  system,  none 
of  these  methods  ensure 

a.  the  stability  of  the  system,  and 

b.  the  dominance  of  the  selected  roots. 

B.    PHILOSOPHY  OF  OPTIMIZATION 

Since  placing  a  pair  of  complex  roots  at  the  desired 
locations  has  no  significance  without  ensuring  their  domi- 
nance and  the  system's  stability,  the  method  proposed  in 
the  first  part  of  the  thesis  provides  these  two  conditions 
by  confining  the  undesired  roots  to  appropriate  region (s) 
of  the  s-plane,  and  then  finds  the  optimum  parameter  values 
to  locate  the  dominant  roots  as  close  to  the  specified 
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1 
locations  as  possible  .   Since  the  undesired  roots  are  not 

placed  at  specific  locations  (they  are  only  confined  to 
stay  in  a  remote  region  of  the  s-plane)  none  of  the  free 
parameters  are  "used  up";   only  their  freedom  is  limited. 
The  problem  is  solved  in  the  "Parameter  Space"  which 
is  a  N  -dimensional  Euclidean  space.  A   performance  index 
(J) ,  which  is  usually  a  measure  of  closeness  of  the  dominant 
roots  to  the  desired  root  locations,  is  defined.   By  mini- 
mizing the  performance  index  with  respect  to  the  selected 
free  parameters  (z)     (see  Section  II BD)    one  obtains  the 
"best  possible  solution"  to  the  design  problem  as  it  is 
specified.   The  design  specifications  yield  some  linear  and 
nonlinear  constraining  equations  and/or  inequalities;   it 
is  assumed  that  these  constraining  equations  define  a  bound- 
ed region  in  the  N  -dimensional  parameter  space.   Some  of 
the  constraining  equations  are  generally  nonlinear  and  it 
may  sometimes  be  impossible  to  find  a  feasible  point  (  a 
point  which  simultaneously  satisfies  all  of  the  constrain- 
ing equations  and  inequalities)  to  minimize  the  performance 
index.   In  practical  terms,  this  means  that  the  design  (or 
compensation)  is  not  always  attainable  by  using  a  certain 


For  the  simplicity  of  explanation  it  is  assumed  here 
that  the  design  is  carried  out  for  a  selected  pair  of 
complex  roots;   but  this  is  not  mandatory.   The  method  is 
generalized  (in  Section  II. B)  and  an  example  is  given 
(Example  3)  for  the  design  of  systems  without  a  pair  of 
selected  dominant  roots. 
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type  of  compensator.   If  this  occurs,  the  designer  should 
change  the  system  configuration  (usually  by  changing  the 
type  of  compensator,  or  by  inserting  additional  compensa- 
tion) . 

C .  CONSTRAINTS 

Performance  specifications  and  available  system  com- 
ponents dictate  the  constraints,  which  can  be  divided  into 
three  general  groups: 

1.  Constraints  due  to  physical  limitations,  such  as 
limitations  on  the  available  system  components. 

2.  Constraints  due  to  static  (or  steady-state)  per- 
formance specifications. 

3.  Constraints  due  to  dynamic  performance  specifica- 
tions. 

The  constraints  in  the  first  two  groups  are  linear  or 
can  be  approximated  by  linear  inequalities.   The  constraints 
in  the  last  group  are  obtained  by  factoring  the  system's 
characteristic  equation  (as  explained  in  the  next  section) 
and  they  are  expressed  as  "n"  linear  and/or  nonlinear  si- 
multaneous algebraic  equations. 

D,  FACTORING  THE  CHARACTERISTIC  EQUATION 

The  denominator  of  the  system's  transmission  function 
(from  input  to  output)  is  called  "the  characteristic  poly- 
nomial", and  when  it  is  equated  to  zero  it  is  called  "the 
characteristic  equation".   The  characteristic  polynomial  of 


15 


th 
the  system  is  a  n   -degree  polynomial  in  the  complex  vari- 
able "s",  and  a  function  of  N   free  system  parameters  (a)  ; 
that  is: 


n 
ch 


,_  (s)  =  2,  ai  sl       *  en. i) 


i=0 


where  a   =1  and  other  coefficients  are  either  known  con- 
n 

stants  or  functions  of  one  or  more  system  parameters. 
Therefore,  the  characteristic  polynomial  can  be  represented 
as 


Pch(s)  =  Pch(^'s)  •  (II'2) 


To  represent  the  desired  and  undesired  roots  explicit- 
ly, it  is  necessary  to  factor  the  characteristic  polynomial 
into  two  parts,   In  this  factored  form  the  characteristic 
polynomial  may  be  written  as 

Pch<S.S>   =  j>d«T.«.s)]  •  LPud'S'S)J  <"-3) 

where  j3  represents  the  coefficients  in  P  ,  (a  polynomial  of 

degree  "n-2"  in  s,  representing  the  undesirable  roots)  and 

P,  is  a  quadratic  in  s  which  represents  the  desired  roots. 

Two  forms  of  P.,  will  be  considered: 
d 

2  2     "• 

P  (<j,lo,s)    =   s   +  2crs  +(<j   +  co),  (If.*4) 
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where   0"   and   to  are  the  real  and  imaginary  parts  of  the 
desired  complex  roots,  and 

P,(C,o>  ,s)  =  s2  +  2Cw  s  +  CO  2,  (II. 5) 

d    n  n      n 

where  C  is  the  damping  ratio  and  CO  is  the  natural  frequen- 
cy of  the  second-order  system. 

The  characteristic  equation  can  be  factored  by  using 
several  methods,  such  as,  by  dividing  the  characteristic 
equation  by  P ,  or  P  ,  and  equating  the  remainder  (in  each 
case)  to  zero.   In  another  method  the  general  forms  of  P, 
and  P  ,  are  multiplied  together  and  the  coefficients  of  the 
product  polynomial  are  equated  to  the  corresponding  coef- 
ficients of  the  characteristic  polynomial;   that  is,  the 
identity  shown  below  is  solved 


Pd  •  Pud  5  Pch  •  (II-6) 


These  processes  yield  n  algebraic  equations,  all  of 
which  must  be  satisfied  to  obtain  the  desired  factored 
form  of  the  characteristic  equation. 

The  left  side  of  Eq.  (II. 3)  contains  N   free  system 
parameters  (a) ,  and  the  right  side  is  a  function  of  n  root 
parameters  (cr,^  ^)  ;   that  is,  there  are  a  total  of  N  +  n 
parameters  available,  but  only  N  of  them  are  independent  — 
the  others  are  related  to  the  independent  ones  through  the 
n  algebraic  constraining  equations  which  are  obtained  during 
the  factoring  process.   Any  N  of  these  parameters  can  be 
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chosen  as  free  variables  (z) ;   the  remaining  n  variables 
are  determined  by  the  constraining  equations  due  to  dynamic 
performance  specifications  (see  Paragraph  II.C(3)  above). 

E.    THE  PERFORMANCE  INDEX  (J) 

When  the  characteristic  equation  is  factored,  that  is, 
after  the  free  variables  (z)  are  chosen  and  the  constrain- 
ing equations  are  found,  a  performance  index  (J)  is  defined 
so  that  when  it  is  minimized  with  respect  to  the  chosen 
variables  (z)  in  the  bounded  region  defined  by  the  con- 
straining equations  and  inequalities,  it  yields  "the  best 
root  locations".   The  best  root  locations  may  mean  that  the 
desired  roots  are  as  close  to  the  specified  locations  as 
possible  (if  the  dominant  roots  concept  is  being  used  in 
the  design) ,  or  to  satisfy  other  criteria  —  as  explained 
below  —  and  that  the  undesired  roots  are  in  remote  regions 
of  the  left-hand  side  of  the  s-plane  (see  Figs.  II. 1  and 
II. 2) . 

The  form  of  the  performance  index  depends  on  the  de- 
sign specifications.   If  the  design  is  to  be  carried  out  by 
locating  a  pair  of  dominant  roots  as  close  to  the  desired 
locations  as  possible,  then  the  performance  index  may  take 
the  form 

J  =  (crd  -  cj)2  +  (0Jd  -  to)2   ,   or  (II. 7) 

J  =  (C.  -  C)2  +  (w   -  a>  )2  (II. 8) 

d  v  n .,    n 

d 

where  the  subscript  "d"  indicates  the  desired  coordinates. 
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Fig.  II. 1.   The  desired  and  undesired  root  location  in 
the  s-plane. 


NOTES : 


Point  (A)  is  an  actual  root  location  at  a 
minimization  step  (J  =  r  2') 

r„  is  a  predetermined  norm  and  it  is  equi- 
valent to  f  in  Fig.  II. 5. 


Fig.  II. 2   The  desired  and  undesired  root  locations  in 
the  s-plane. 

(Point  (A)  is  an  actual  root  location  at  a 
minimization  step) . 
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Sometimes  design  specifications  allow  the  designer  to 
use  a  simpler  performance  index.   For  example,  only  the 
first  terms  of  Eq.  (II. .7)  or  (II. 8)  may  be  used.   For  some 
purposes  extra  terms  may  be  added  to  the  performance  index. 
For  example,  it  may  be  desired  to  confine  the  desired  roots 
to  two  small  regions  around  the  desired  root  locations  and 
the  real  part  of  the  closest  undesired  root  (to  the  origin) 
is  maximized  by  minimizing  the  performance  index. 

J  =  -q   ,   (q  >  0)  (II. 9) 

where  q  is  the  magnitude  of  the  real  part  of  the  closest 
undesired  root.   The  negative  sign  in  (II. 9)  converts  the 
maximization  problem  into  a  minimization  problem. 

•  For  different  kinds  of  dynamic  performance  specifica- 
tions, different  performance  indices  can  be  defined  and 
minimized  with  respect  to  the  chosen  free  parameters.   If 
the  dynamic  performance  specifications  are  very  tight  (that 
is  certain  rise  time,  overshoot,  settling  time,  frequency 
of  the  transients,  etc.,  are  desired)  it  is  best  to  use  the 
dominant-roots  concept  (Eq.(II.7)  or  (II. 8))  for  the  per- 
formance index. 

F.    PROGRAMMING 

Although  it  is  possible  to  solve  simple  problems  by 
hand  calculations,  in  general,  a  computer  program  capable 
of  minimizing  a  nonlinear  function  with  linear  and  nonlinear 
constraints  is  necessary  for  higher-order  systems.   For  the 
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examples  considered  in  this  part  of  the  thesis  SUBROUTINE 
BOXPLX  (see  Refs.  1  and  4)  has  been  used  in  the  minimiza- 
tion with  good  results.   A  computer  program  which  is  written 
in  FORTRAN  IV,  and  used  for  the  minimization  part  of  Ex- 
ample -  4  is  given  at  the  end  of  the  thesis  (see  COMPUTER 
PROGRAM  I)  ., 

G.    EXAMPLE  1 

A  linear  third-order  feedback  control  system  is  shown 
in  block  diagram  form  in  Fig.  11.3(a).   A  step-by-step  pro- 
cedure to  compensate  the  system  for  the  desired  transient 
response  by  using  the  method  "Optimization  for  the  Best 
Root  Locations  in  the  s-Domain"  is  given  below. 

1.  Interpret  the  Desired  Dynamic  Performance  Specifi- 
cations 

The  desired  dynamic  performance  specifications  are 
usually  given  as  maximum  overshoot,  rise  time,  settling  time, 
frequency  of  the  transients,  etc.,  and  it  is  usually  best 
to  find  a  pair  of  complex  root  locations  in  the  s-domain  to 
represent  these  specifications.   For  this  problem  it  is  as- 

tat 

sumed  that  a  pair  of  cqmplex  roots  located  at  (J-.   =  0.5  and 
U),  =  ±0.5  would  yield  the  desired  dynamic  response. 

2.  Check  to  Determine  if  the  Desired  Roots  are  on  the 
Root  Locus  of  the  Uncompensated  System 

The  designer  gains  considerable  insight  into  the 
design  problem  by  just  sketching  the  root  locus;   this  also 
helps  in  deciding  on  the  type  of  compensation  that  should 
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be  used,  if  different  choices  are  available.   The  root 
locus  for  the  uncompensated  system  is  shown  in  Fig.  11.3(b), 
and  the  desired  locations  are  not  on  the  root  locus. 

3 .  Decide  on  the  Type  of  Compensation 

A  cascade  compensator  will  be  used.   The  system 
with  the  cascade  compensator  is  shown  in  Fig.  11.4(a). 

4.  Find  the  Transmission  Function  and  the  Character- 
istic Equation 

The  compensated  system  becomes  fourth  order  (n=4), 
and  there  are  three  free  system  parameters  (N  =3)  —  the 
pole  location  (p) ,  the  compensation  ratio  (k)  of  the  com- 
pensator, and  the  system  gain  (K) .   The  vector  representing 
the  free  system  parameters  is 

a   =  (k  P  k)t  .  (n.io) 

The    transmission    function   of    the   compensated   system   is 

T(s)   =  |isl  .   (Vk)  (a  ±  kp) ,        (II.n) 

R{s)         s(s+l)     (s+p)    +    (K/k) (s+kp) 


and  the  characteristic  equation  is 

s4  +  (2+p)s3  +  (l+2p)s2  +  (p+K/k)s  +  Kp  =  0  .  (11.12) 


5.   Factor  the  Characteristic  Equation 

Using  one  of  the  methods  described  in  Section  II. D 
factor  the  characteristic  equation  and  find: 

a.   The  (n-2)    =  2   -degree  polynomial  (P  ,)  to 
represent  the  undesired  roots. 
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J*> 


(a)  System 


Desired  root  locatio 


#- 


s=-l 


S=-crJ=.o5 


(b)   Uncompensated  system's  root  locus 


Fig.  II. 3.   EXAMPLE  1 


*-% 


(1/k) (s+kp) 


s(s+l)g 


DESIRED  ROOT 
C(s)   LOCATION 


(a)  Compensated  system 


S-.94 

(b)  Root  locus  for  the  compensated  system 

Fig.  I I. 4.  EXAMPLE  1 


Sj-O.S 


*/0 


* 
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Pud(j?,s)  =  s2  +  2as  +  b  (11.13) 

m 

where  £  =  (a  b) 

b.  A  quadratic  in  s  (P  )  to  represent  the  desired 
roots  (the  roots  which  will  be  made  dominant  by  optimiza- 
tion) . 

P  (cr,a>,s)  =  s   +  2d    s  +  (CT2  +  co2)  (11.14) 

c.  The  type  -3  constraining  equations  by  using 
Eq.  (II. 6) 

Pch(~'s)  =  [Pd^^'S)  J  '  LPud(J?'s)]  *  (II- 6) 

2     ? 

For  notational  convenience  let  c  =  (cr   +  a)  )  , 


s4  +  (2+p)s3  +  (l+2p)s2  +  (p+K/k)s  +  Kp  = 


(11.15) 


4  3  2 

s   +  2(a+cr)s   +  (c+4acr+b)  s  +  2  (ac+bcr)  s  +  be  . 


Equate  the  coefficients  of  the  polynomials  of  the  above 
identity  and  obtain  the  n(=4)  constraining  equations 

p  +  K/k  =  2(ac  +  bd)  (11.16) 

Kp  =  be  (11.17) 

1  +  2p  =  4acr  +  b  +  c  (11.18) 

2  +  p  =  2(a  +  cr)  (11.19) 


There  are  total  of  N   +  n  (=  7)  variables  (p,K,k,a,b, 
o^O))  ,  and  any  N   (=  3)  of  them  can  be  chosen  as  the  free 
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variables  to  constitute  the  z-vector  of  the  minimization; 
the  remaining  n  (=  4)  will  be  related  to  the  free  variables 
by  the  constraining  equations.   After  selecting  the  free 
variables,  the  constraining  equations  should  be  rearranged 
so  that: 

(1)  the  left  side  of  each  equation  is  one  of 
the  dependent  variables,  and 

(2)  the  right  side  of  each  equation  must  con- 
tain only  the  free  variables  and  the  dependent  variables 
defined  in  the  preceding  equations. 

To  avoid  substitution,  (7,  a,  and  c  are  chosen  as  the 
free  variables,  hence 

z  =  (a  cr  c)T  ,  .     (11.20) 

and  the  constraining  equations  are  solved  for  the  remaining 
four  variables 

p  =  2  (a  +  <J  -  1)  (11.21) 

b  =  l  +  2p-c-  4a<7  (11.22) 

K  =  bc/p  (11.23) 

k  =  K/(2  (ac  +  b(T)  -  p)  ,                      (11.24) 

6.   State  the  Additional  Constraints 

a.   The  compensator  will  contain  only  passive  ele- 
ments;  hence, 

p  >  0  (11.25) 
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b.  7A  physical  constraint  is  assumed  on  the  sys- 
tem gain 

0.4  ^  K  £  0.5  (11.26) 

c.  The  compensation  ratio  of  a  cascade  compensa- 
tor should  not  be  greater  than  ten. 

0.1  £  k  ^  10  (11.27) 

d.  To  ensure  complex  dominant  roots, 

a2  -  c  <  0  (11.28) 

e.  To  confine  the  undesired  roots  to  the  left  of 
s  =  -(T,   =  -1.5  line, 

a  -  (a2  -  b)^  *  a    =  1.5  (11.29) 

f.  All  free  and  dependent  parameters  must  be  pos- 
itive;  hence 

a   >   0 

c  >  0  (11.30) 

b  >  0 

7 .   Define  a  Performance  Index 

For  this  problem  the  performance  index  defined  by 
Eq.  (II. 7)  is  suitable, 

J  =  (<Jd  -  cr)2  +  (u>d  -  w)2  (II. 7) 
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8.  Minimize  the  Performance  Index 

Find  the  best  root  locations  and  the  corresponding 

optimum  system  parameters  (a*)  by  minimizing  the  performance 

T 
index  with  respect  to  the  chosen  free  parameters  z  =  (a  cr  c) 

in  the  bounded  parameter  space  defined  by  the  constraining 

inequalities  (11.25)  -  (11.30)  and  the  constraining  Eqs. 

(11.21)  -  (11.24). 

9.  Relax  the  Necessary  Boundaries  if  Possible 

If  at  the  end  of  a  minimization  process  the  select- 
ed roots  cannot  be  located  in  a  close  neighborhood  of  the 
desired  locations  (determined  by  a  preselected  norm  —  see 
Eq.  (11.31)  below),  and  the  result  is  a  constrained  minimum 
(that  is,  at  least  one  of  the  variables  is  on  the  limiting 
boundary)  then  the  limiting  boundary (s)  should  be  relaxed 
(if  possible)  and  the  minimization  should  be  repeated.   In 
this  example  the  a.   value  was  the  only  boundary  which 
limited  the  minimization  process.   The  boundary  on  the  un- 
desired  roots  was  relaxed  toward  the  origin  by  0.1  steps 
until 

_l/2  *  €  =  10"4  (11.31) 

J 

was  obtained.   After  each  relaxation,  SUBROUTINE  BOXPLX  was 
called  and  the  performance  index  was  minimized  again  to  ob- 
tain optimum  parameter  values  for  the  new  boundaries.   The 
results  obtained  at  the  end  of  each  minimization  are  tabu- 
lated in  Table  II. 1.   Relaxation  was  stopped  when  the 
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undesired  root  (#1)  was  no  longer  on  the  or..   boundary  (see 
last  line  of  Table  II .1);   at  the  same  step  the  stopping 
criterion  on  the  performance  index  (Eq.  (11.31))  was  also 
satisfied. 

The  repeated  minimization  process  and  boundary  relax- 
ation were  achieved  automatically  in  the  main  program  (see 
the  computational  flow  chart  —  Fig.  II. 5). 

10.  Interpret  the  Results 

At  the  first  minimization  step  the  undesired  roots 
were  kept  to  the  left  of  s  =  -<T,  .  =  -1.5  line.   The  desired 
roots  were  probably  dominant  (since  the  transient  due  to 
the  undesired  real  root  settles  down  in  1/6    of  the  time 
of  the  transients  caused  by  complex  roots) ,  but  the  domi- 
nant roots  were  not  close  to  the  desired  locations  and  the 
transient  response  would  not  be  satisfactory.   At  the  last 
step  of  the  relaxation  the  desired  roots  were  exactly  locat- 
ed but  their  dominance  was  questionable  since  an  undesired 
real  root  was  located  close  to  the  desired  roots. 

To  make  the  final  decision  on  the  parameter  values, 
residues  can  be  calculated;   this  would  certainly  remove 
the  doubts  about  the  dominance  of  the  desired  roots,  but 
since  a  digital  computer  is  assumed  to  be  available  and 
digital  simulation  of  the  system  avoids  the  lengthy  residue 
calculations,  it  is  best  to  obtain  the  time  response  of 
the  system  and  compare  it  with  the  desired  second-order 
response. 
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LIMITING 
VALUE    FOR 
THE    UNDE- 
SIRED 
ROOTS  (-CI J 

UNDESI 
RC 

#1 

RED 
)OTS 

#2 

DESIR 
R 

Real 
Part 

ED 
OOTS 

Imag. 
Part 

SYSTEI^ 

k 

1   PA RAM 
P 

ETERS 
K 

J2 

-1.5 

-1.5 

-27.15 

-0.249 

±0.496 

0.24 

27.15 

0.46 

0.25 

-1.4 

-1.4 

-27.56 

-0.297 

±0.469 

0.1 

27.55 

0.43 

0.205 

-1.3 

-1.3 

-    9.85 

-0.336 

±0.492 

0.17 

9.83 

0.463 

0.165 

-1.2 

-1.2 

-13.21 

-0.386 

±0.494 

0.1 

13.18 

0.473 

0.114 

-1.1 

-1.1 

-11.23 

-0.429 

±0.499 

0.1 

11.19 

0.479 

0.07 

-1.0 

-1.0 

-   9.733 

-0.473 

±0.497 

0.1 

9.68 

0.474 

0.02 

-0.9 

-0.94 

-  9.525 

-0.5 

±0.5 

0.1 

9.467 

0.474 

0.00008 

TABLE    II.l 
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INPUT:  Upper  and  lower  limits  on  the  free  variables, the  starting  point 


SUBROUTINE  BOXPLX:  Finds  the 
minimum  in  the  given  bounded 
region, by  using  COMPLEX  method 


Continue  with 
the  normal 
minimization 
process 


zc 


n 


'Generate'^  A — 
'another  ■'  / 
Ipoint 


£ke=o 
1  v_ 


T 


FUNCTION  KE(X) 
Checks  for  all  constraints. 
If  a  point  is  in  the  de-   _^. 
signated  region  (feasible 
point),  set  KE  =  0.  Other- 
wise  set  KE  =  -1. 


YES 


FUNCTION  FE(X) 
Calculates  the  value  of 
the  function  to  be  minim- 
ized (value  of  J) 


1 


A  MINIMUM  IS  FOUND 

PRINT:  Of*,  and  the  values  of 

the  other  necessary  variables 


Are  the  selected  roots  in  a 
close  neighborhood  of  the 
.selected  locations  ? 


"YES 


any  of  the  va- 
riables on  the 
preset 
boundaries? 


YES 


O 


IF  POSSIBLE 
RELAX  THE 
BOUNDARY (S) 
Limiting  the 
minimization 
process  in 
KE(X) 

Are  all 
possible  relaxa- 
tions completed? 


J 


PRINT:  "THE  BEST  ROOT  LOCATIONS 
ARE  FOUND" 

and 
Of*  ,  other  necessary  informa- 
tion  


c 


the  starting  point  to 
zn .        To  start  a  new 
Minimization . 


SET: 


STOP 


} 


6 


€    is   the  predetermined  positive  constant    (Stopping  criterion) 

Fig.    II. 5.    Computational    flow   chart    for    the  method    "Opti 
mization    for    the   Best   Root   Locations    in    the 
s -Domain" . 
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Fig. II. 7.  EXAMPLE  1. 


System's  time  response  (A)  with 
optimum  parameter  values  and 
desired  response  (B)  to  a  unit 
step  input. 
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Fig.  11^6  is  the  response  (to  a  unit  step  input)  for 

the  first  line  of  Table  II. 1;   the  desired  second-order 

response  is  also  shown.   Fig.  II. 7  shows  the  time  response 

of  the  system  (to  a  unit  step  input)  for  the  parameter 

values  tabulated  as  the  last  line  of  Table  II. 1,  where  the 

complex  roots  were  exactly  located  at  the  desired  locations 

(CT  =  0",  =  0.5  and  to  =  to ,  =  ±  0.5)  but  a  real  root  was  close 
d  d 

to  the  desired  roots  (at  s  =  -0.94).   Since  the  compensated 
system's  response  is  very  close  to  the  desired  response, 
the  optimum  parameter  values  tabulated  as  the  last  line  of 
Table  II. 1  yield  a  very  satisfactory  design. 

H.    EXAMPLE  2 

The  linear  fourth-order  system  shown  in  Fig.  118(a)  is 
given.   Steady-state  performance  specifications  require  the 
system  gain  to  be  equal  to 
R(s) 


K 

s(s+l) (s+5) (s+10) 


C(s)  R(s)   + 


+ 


K 

s(s+l) (s+5) (s+10) 


C(s) 


K  s  +  K.  s 
a      t 


(a)  The  uncompensated  system     (b)  The  compensated  system 

Fig.  I I. 8.  EXAMPLE  2 

(or  greater  than)  1000  (K  ^  1000) ,  but  the  root-locus 
study  shows  that  the  system  is  unstable  if  K  >  193.23. 
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1.   The  dynamic  performance  specifications   can  be  met 

by  locating  a  pair  of  complex  (dominant)  roots  at  (J,  =  1   ■ 

and  «,  =  ±  2. 
d 

2„   The  desired  root  locations  are  not  on  the  uncom- 
pensated system's  root  locus,  and  the  system  gain  is  much 
less  than  the  desired  value  (for  the  part  of  the  root  locus 
in  the  left  hand  side  of  the  s-plane) ;   therefore,  some  kind 
of  compensation  is  needed. 

3.  A  combination  of  tachometer  and  acceleration  feed- 
back will  be  used  to  compensate  the  system.  The  configura- 
tion of  the  compensated  system  is  shown  in  Fig.  118(b). 

4.  The  transmission  function  of  the  compensated  system 
is 

T(S)  =  R(ff  =  ~4 ~3 ~ 2 '  (II-32) 

K    '  s  +  16s  +(65+KK  )s  +(50+KK.)s+K 

a7    x      t 

and  the  characteristic  equation  is 

s4  +  16s3+  (65+KK  )s2+  (50+KK.)s  +  K  =  0  .       (11.33) 

a  c 

5.  The    desired-rooLo  quadratic    is 

Pd(cr,   CO,    s)    =   s2    +   2crs    +   c  (11.34) 


where 


^2  2 

c   =  a     +  <jc 


Paragraph  numbers  indicate  the  corresponding  steps  in 
EXAMPLE  1. 
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ud 
where 


The  undesired-roots  polynomial  is 

P  ,(£,  s)  =  s2  +  2as  +  b  (11.35) 


/9  -  (a  b)T  . 


Eq.  (II. 6)  is  used  to  factor  the  system's  character- 
istic equation, 


s4+  16s3+  (65+KK  )s2+  (50+KKjs  +  K  = 

a  t 

4  3  2 

s  +  2(a+cr)s  +  (b+c+4acr)  s  +  2(bc+ac)s  +  be. 


(11.36) 


By  equating  the  coefficients  of  the  corresponding 
terms  of  identity  (11.36), 

a   +   (J   =   8 

b   +   c    +   4aCT   =   65    +   KK 

a  (11.37) 

2(bcr  +   ac)    =   50   +  KK 

be    =   K 

The  z-vector,  representing  the  independent  variables 
in  the  minimization  process,  is  chosen  to  be 

z  =  (a  c  K)T  .  (11.38) 

The  constraining  equations  are 
a  =  8  -<7 
b  =  K/c 


K  =    (b   +   c    +   4a<7    -    65)/K 
K   =    (2     (bCT    +    ac)     -    50)/K 
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(11.39) 


6.   All  the  variables  -must  be  positive  (or  zero),  and 
the  other  constraints  are 


K  £  1000 


0  ^  K   £  10 
a 

0  ^  K.  ^  10 
t 


a   -  c  <  o 


2     -, 


a  -  (a  -b  )   *  crn  ,  =  6 


7.   The  performance  index  is 


J  =  (1,-CT)2  +  (2. -CO)2 


(11.40) 


(II. 40a) 


(11.41) 


8.  SUBROUTINE  BOXPLX  was  again  used  for  the  minimiza- 
tion. 

9.  None  of  the  variables  was  on  the  limiting  bound- 
aries;  therefore,  no  relaxation  was  necessary.   The  result 
of  the  minimization  is  tabulated  in  Table  II. 2. 

10.   The  compensated  system's  response  (to  a  unit  step 
input)  and  the  desired  response  are  shown  in  Fig.  II. 9. 


LIMITING 

UNDESIRED  ROOTS 

DESIRED  ROOTS 

SYSTEM 

■ 

'  VALUE 

PARAMETERS 

FOR  THE 

UNDE- 

SIRED 

ROOTS 

Real 

Imag, 

Real 

Imag. 

<-<Tlt) 

part 

part 

part 

part 

K* 

K*    1 
a 

K* 

-6 

-6.999 

-199.585 

-1.001 

±2.002 

1000 

.1676  1 
i 

.4197 

TABLE  I I. 2 
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(a)  ;  //J  -(b) 


1. 


2. 


TIME 


5* 


Fig. II. 9.  EXAMPLE  2. 


Time  response  (A)  of  optimally  com- 
pensated system  and  the  desired  re- 
sponse (B)  (to  a  unit  step  input). 
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I.    EXAMPLE  3 

The  linear  system  given  in  Example  2  (and  shown  in 
Fig,  11.8(a))  will  be  compensated  for  the  fastest  response 
(the  shortest  rise  time  and  the  shortest  settling  time). 
The  steady-state  requirements  and  the  type  of  compensation 
to  be  used  are  exactly  the  same  as  in  Example  2. 

Since  the  only  dynamic  system  specification  is  "obtain 
the  fastest  response",  the  actual  form  of  the  transient  re- 
sponse and  the  location  of  the  dominating  roots  (whether 
they  exist  or  not)  are  of  no  interest  to  the  designer.   This 
example  shows  that  the  dominant-root  idea,  although  very 
useful  for  many  design  problems,  is  not  mandatory  for  the 
system   design  procedures. 

As  mentioned  earlier  (in  Section  II. A)  the  s-domain 
design  methods  try  to  relocate  the  roots  of  the  system's 
characteristic  equation  to  obtain  a  suitable  transient  re- 
sponse.  Most  of  these  methods  ignore  the  magnitude  effects 
(to  avoid  the  residue  calculations),  and  they  try  to  clas- 
sify the  transient  response  as  "the  fast-varying  part",  and 
"the  slowly-varying  part".   This  procedure  depends  on  the 
assumption  that,  "magnitude-wise  all  of  the  roots  have  the 
same  importance  in  the  transient  response". 

With  the  same  assumption  one  may  try  to  obtain  "the 
fastest  response"  by  making  the  real  parts  of  all  the 
roots  as  large  as  possible,  because  the  rise  time  and  the 
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settling  time  of  a  transient  response  caused  by  a  root  is 
inversely  proportional  to  the  real  part  of  that  root. 

The  solution  of  the  problem  contains  exactly  the  same 
steps  as  the  previous  example  (see  Eqs.  (11.32)  -  (11.41)). 
Only  the  last  two  equations  should  be  changed,  specifically 
Eqs.  (II. 40a)  and  (11.41). 

Eq.  (II. 40a)  was  a  constraining  equation  which  ensured 
the  confinement  of  the  undesired  roots  to  the  left  of  the 
s  =  -(J,   line  for  the  case  of  two  real  undesired  roots.   In 
this  problem  it  takes  the  form 


2 


L 


a  -  (a^  -  b)  2  ^  cr  .  (11.42) 

If  the  undesired  roots  are  complex  ,  the  equation 


a  ^  a  (11.43) 

serves  the  same  purpose. 

These  two  constraining  equations  ensure  that  the  magni' 
tudes  of  the  real  parts  of  the  undesired  roots  are  greater 
than  the  magnitudes  of  the  real  parts  of  the  desired  roots. 
The  performance  index  takes  the  simple  form 

J  =  -cr   ,   (or  *  0)  ,  (11.44) 

because  it  is  only  necessary  to  maximize  the  real  parts  of 
the  desired  roots. 


For  this  problem  all  roots  are  undesirable,  but  just 
for  purposes  of  terminology  the  names  used  in  Example  2 
are   retained.   The  roots  of  the  polynomial  PU(j  are  still 
called  "the  undesired  roots"  and  the  roots  of  P^  are  cal- 
led "the  desired  roots". 
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The  results  of  the  minimization  are  shown  in  Table 
II „ 3,  and  the  time  response  to  the  unit  step  input  is  shown 
in  Fig.  11.10.   The  desired  response  for  Example  2  is  also 
plotted  in  the  same  figure  for  easy  comparison. 


UNDESIRJ 

Real 
part 

3D  ROOTS 

Imag. 
part 

DESIRED  ROOTS 

Real     Imag. 
part    part 

SYSTEM  PARAMETERS 

v  '     k*      K* 
a       t 

-4. 

±34.8 

-4. 

±3.57 

1000. 

.0625 

.458 

TABLE  II. 3 

J.    EXAMPLE  4 

In  some  design  problems  it  may  be  necessary  to  insert 
a  lag  compensator  which  has  a  dipole  (a  pole  and  a  zero) 
very  close  to  the  origin.   In  this  case  some  undesired 
root(s)  may  be  located  closer  to  the  origin  than  the  de- 
sired roots.   Problems  of  this  type  may  be  solved  by  using 
the  normal  step-by-step  procedure  given  in  Example  1;   if 
this  is  done,  an  acceptable  answer  may  be  obtained  by  re- 
laxing the  limiting  value  for  the  undesired  roots  (on , )  by 
a  considerable   amount  toward  the  origin  (as  in  step  9  of 
Example  1) .   If  the  designer  can  gain  enough  insight  to 
the  problem  (probably  at  step  2  of  the  normal  step-by-step 
procedure)  he  may  directly  select  OV   =  0  (to  ensure  the 
stability  of  the  system)  and  save  some  computer  time  by 
avoiding  the  relaxation  steps  for  cr,  .  . 
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Fig.  11.10.  EXAMPLE  3.   The  fastest  time  response  (A)  and 

the  desired  response  (B)  for 
EXAMPLE  2. 
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In  another  approach  the  designer  may  identify  and  con- 
fine some  of  the  undesired  roots  to  a  close  vicinity  of  the 
origin  while  keeping  the  other  undesired  roots  in  a  remote 
region  of  the  s-plane.   In  the  design  example  given  below 
this  method  is  used  just  to  explain  a  variation  from  the 
main  method  given  in  Example  1. 

A  third-order  system  (shown  in  Fig.  11.11(a))  is  given. 
Steady-state  performance  requirements  specify  that  the  sys- 
tem gain  (K)  must  be  greater  than  425  (the  stability  limit 
of  the  uncompensated  system) .   The  step-by-step  solution  is 
given  below;   paragraph  numbers  indicate  the  corresponding 
steps  in  Example  1. 

R 


Rfe) 

K 

C 

:(s) 

-^y 

\ 

s(s+6)2 

— i 

i — »• 

*%•: 


1      s+kp 
k      s+p 

K 

c( 

s) 

> 

s(s+6)2 

(a)  Uncompensated  system 


(b)  Compensated  system 


Fig. II. 11.  EXAMPLE  4 


1.  Dynamic  performance  specifications  require  a  pair 
of  dominant  roots  located  at  cr,  =  1,  and  to..  =  ±  2. 

2.  -  3.   The  designer  decides  on  a  cascade  lag  com- 
pensator by  sketching  the  root  locus  of  the  uncompensated 
system.   The  compensated  system  is  shown  in  Fig.  II. 11(b). 
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4.      The    transmission    function   of    the  compensated   system 
is 


T(s)    =  (K/k)  (s+kp) f        (11.45) 

s4+(12+p)s3+(36+12p)s2+(36p+K/k)s+Kp 


and    the   characteristic    equation   is 

Pch(tt,  s)=  s4+CL2+p)s3+(36+12p)s2+(36p+K/k)s+Kp   =   0 

(11.46) 
where 

a  =    (k  p  K)T   . 

5.   The  quadratic  for  the  desired  root  is 

Pd((7fU5,s)  =  s2  +  2(7  s  +  c  (11.47) 

where 

a   2     2 
c    =   G      +  lc      . 

For  this  prcblem  the  characteristic  equation  of  the 

system  cannot  have  more  than  two  complex  roots;   since  the 

designer  is  interested  in  the  part  of  the  root  locus  where 

the  desired  roots  are  complex,  then  the  undesired  roots 

must  be  real.   For  this  reason  the  polynomial  representing 

the  undesired  roots  can  be  written  as 

P  d(£,s)  =  (s+a)-(s+b)  =  s2+  (a+b)s  +  ab        (11.48) 

The  characteristic  equation  is  factored  by  using  iden- 
tity (II. 6) 

P  ,  =   P ,  P  ,  (II. 6) 

ch     d   ud 
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s4+(12+p)s3+    (36+12p)s2+    (36p+K/k)s    +  Kp    =  (11.49) 

4  3  2 

s    +(a+b+2(J)s    +  (ab+c+2a(J+2bcr)  s    +  (2abcr+ac+bc)  s    +   abc, 


and  (J,    p,    and  b   are   chosen   as    the  variables    in    the  minimiza- 
tion process;    in   other  words, 

z  =  (a  p  b)T  (ii.50) 

The  constraining  equations  can  be  written  (from  Eq. 
(11.49))  as 

a  =  12+p-b-2(J 

c  =  36  +  12p  -  2a(a+b)  -  ab 

(11.51) 
K  =  abc/p 


k  =  K/(-36p  +  2aba  +  c(a+b)) 


6.  All  of  the  variables  must  be  greater  than  or 
equal  to  zero,  and  the  other  constraining  inequalities  are 

K  ^  425 

1  <  k  ^  10  (11.52) 

b  >  0",  ,_  ,    limit  for  the  remote  undesired  root 
lt-1 

0  <  a  ^  O.^  0  limits  for  the  close  undesired  root. 

7.  The  performance  index 

J  =  (l.-CT)2  +  (2.-W)2  (11.53) 

is  used  for  minimization  by  SUBROUTINE  BOXPLX. 
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9.   Some  suitable  logic  to  relax  the  boundaries  was 
inserted  in  the  main  program.   The  remote  undesired  root 
(at  s  =  -b)  is  confined  to  the  left  of  the  s  =  -G        _  line 
and  the  close  undesired  root  (at  s  =  -a)  is  confined  be- 
tween the  s  =  0  and  s  =  -0",  .  ~  lines,  where  cr,  ,  ,  =8,  and 

lt-2  lt-1 

(T,  ,_2  =  1.5.   The  first  limiting  line  would  be  shifted  to 
the  right  and  the  second  one  shifted  to  the  left,  but  no 
relaxation  is  done  because  the  minimum  is  not  on  the  enclos- 
ing boundaries.   The  result  of  the  minimization  is  tabulated 
as  Table  II. 4. 

10.   The  compensated  system's  response  (to  a  unit  step 
input)  and  the  desired  response  are  shown  in  Fig.  11.12. 
Although  the  desired  roots  are  located  exactly  at  the  selec- 
ted locations,  they  are  obviously  not  dominant  becasue  the 
two  time  responses  are  not  close  to  each  other. 


LIMITIN 
FOR  THE 
SIRED  R 

-ait-l 

G  VALUES 

UNDE- 
OOTS 

"ait-2 

UNDESIRED  ROOTS 
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DESIRE 

Real 
part 

D  ROOTS 

Imag. 
part 

SYS' 
PAR. 

K 

rEM 

^MEl 
P 

'ERS 

-8. 

-1.5 

-1.38  j  -8.723 

1. 

±2. 

563.  j. 107 

1 

7.952 

TABLE  I I. 4 

K.    CONCLUDING  REMARKS 

The  method  "Optimization  for  the  Best  Root  Locations  In 
the  s-Domain"  provides  fast,  direct  and  efficient  solution 
to  the  linear  system  design  problems  of  classical  control. 
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Fig.  11.12.  EXAMPLE  4.   Time  response  (A)  of  the  compensated 

system  with  a  lag  compensator  and 
the  desired  response  (B)  to  a  unit 
step  input. 
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Systems  with  more  than  two  free   parameters  can  be  handled 
easily. 

For  high-order  systems,  algebraic  manipulations  neces- 
sary to  factor  the  characteristic  equation  may  get  too  in- 
volved.  If  the  number  of  free  parameters  is  very  small 
compared  to  the  order  of  the  system,  these  algebraic  mani- 
pulations may  become  prohibitively  time  consuming. 
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III.   OPTIMIZATION  FOR  THE  BEST  RESPONSE  IN  THE 
TIME-DOMAIN  —  APPLICATIONS  TO  LINEAR 
SYSTEM  DESIGN  PROBLEMS 

A .    GENERAL 

Most  of  the  methods  for  the  design  of  dynamic  systems 
are  inherently  indirect,  s-domain  methods,  for  example, 
try  to  adjust  the  root  locations  of  the  system  character- 
istic equation  to  obtain  a  suitable  transient  response. 
For  a  high-order  system  it  may  be  very  difficult  to  decide 
on  a  set  of  suitable  root  locations  which  defines  a  desir- 
able transient  response;   for  this  reason,  the  final  step 
of  an  s-domain  design  procedure  is  usually  the  simulation 
of  the  system,  to  check   if  the  new  root  locations  yield 
an  acceptable  transient  response. 

The  method  presented  in  this  chapter  eliminates  the 
intermediate  steps  and  directly  tries  to  shape  the  tran- 
sient response  of  the  system.  This  method  is  called  "Op- 
timization for  the  Best  Response  in  the  Time-Domain". 

The  first  step  of  the  method  is  selecting  a  "MODEL 
RESPONSE"   by  interpreting  the  dynamic  performance  speci- 
fications.  Then  the  state  equations  of  the  system  are 
integrated  with  respect  to  time  with  a  given  input 


1 

The  word  "response"  implies  "time  response"  even  when  it 
is  used  alone. 
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function,  and  the  output  of  the  system  is  obtained  and 
sampled  at  suitable  time  intervals.   These  samples  are  com- 
pared with  the  desired  values  of  the  response  (defined  by 
the  model  response) .   A  performance  index,  which  measures 
the  deviation  of  the  sample  values  from  the  desired  values 
(from  zero  time  to  the  expected  end  of  the  settling  period 
of  the  system)  is  defined.   Minimization  of  this  performance 
index  with  respect  to  the  free  system  parameters  yields  the 
best  time  response  of  the  system,  for  the  given  specifica- 
tions and  the  system  configuration. 

In  this  chapter,  the  simple  philosophy  of  the  method 
is  described  and  the  method  is  applied  to  the  linear  sys- 
tem design  problems  which  were  solved  in  Chapter  II  by 
using  the  method  "Optimization  for  the  Best  Root  Locations 
in  the  s-Domain". 

B.    SELECTING  THE  MODEL  RESPONSE 

Selection  of  the  model  response  is  determined  by  the 
dynamic  performance  specifications.   If  the  system  is  to  be 
operated  always  with  a  specific  input  signal,  then  this 
time  function  (input)  can  also  be  used  as   the  model  re- 
sponse for  the  comparison,   A  step  input  (test  signal)  is 
very  common  in  dynamic  system  design  procedures.   The  same 
step  function  may  also  be  used  as  the  model  response;   if 
this  is  done,  however,  it  may  be  necessary  to  insert  some 
constraints  on  the  states  of  the  system  to  fulfill  other 
dynamic  performance  specifications,  such  as  peak  overshoot, 
etc.  (see  Examples  7  and  9). 
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If  the  dynamic  performance  specifications  are  too 
tight  (that  is,  if  expected  rise  time,  peak  overshoot,  set- 
tling time,  frequency  of  oscillations  are  given  within  some 
small  tolerance  limits)  a  second-order  model  system  may  be 
used  to  represent  the  desired  model  response  easily  and 
efficiently.   The  output  of  the  model  system  is  obtained 
by  integrating  the  state  equations  of  the  model.   Using  a 
second-order  system's  output  as  the  model  response  is,  of 
course,  equivalent  to  the  dominant-roots  concept  in  the 
s-domain. 

C.    CONSTRAINTS 

The  constraints  on  the  free  system  parameters  are 
usually  linear  or  can  be  approximated  by  linear  inequali- 
ties --  this  simplifies  the  computer  programming  part  of 
the  design  procedure. 

There  may  be  some  constraints  on  the  state  of  the  sys- 
tem due  to  inherent  nonlinear  characteristics  of  the  system 
components.   These  constraints  can  be  represented  by  linear 
or  nonlinear  inequalities  and  they  do  not  introduce  much 
difficulty  to  the  design  procedure  (see  Chapter  IV) . 

To  obtain  some  design  characteristics,  additional  con- 
straints on  the  states  of  the  system  may  be  introduced  by 
the  designer.   For  example,  in  one  design  procedure  the_ 
peak  overshoot  of  the  system  can  be  kept  below  a  certain 
level  by  inserting  suitable  constraints  on  the  output 
state.   In  this  thesis,  penalty  functions  are  employed  to 
handle  this  type  of  constraint  (see  Examples  7  and  9) . 
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D.    DEVELOPMENT  OF  THE  METHOD 

The  state  equations  of  an  n   -order,  single-input  sys- 
tem can  be  written  in  the  form 


x(t)  =  Ax(t)  +  br(t) 


(III.D 


where  x(t)  is  the  state  n-vector  and  the  "dot"  over  a  vari- 
able is  used  to  indicate  the  time  derivative  (d/dt) .   A  is 
an  n  x  n  matrix,  r(t)  is  the  scalar  time  function  represent- 
ing the  input  to  the  system  and  b  is  an  n-vector. 
The  output  equation  is 


c  (t)  =  a   x(t)  +  h   r(t) 
s      ~   ~       o 


(III. 2) 


where  c  (t)  is  the  output  of  the  system,  a  is  an  n-vector 

and  h   is  a  scalar, 
o 

Assume  that  the  design  is  to  be  carried  out  by  using 
a  second-order  model  response.   The  state  equations  for  the 
second-order  model  system  can  be  written  as 


y(t)  = 


-to 


n 


a 


d  n 


y(t)  + 


10 


0 

2 
n 


r(t)    (III. 3) 


This  step  may  take  different  forms  depending  on  "the  choice 
of  the  model  response.   If  the  input  time  function  (r(t)) 
is  to  be  used  as  the  model  response,  for  example,  no  extra 
effort  is  necessary  to  generate  the  model  response. 
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where  y(t)  is  the  state  vector  of  the  model  system,  C,  re- 

rv/  a 

presents  the  desired  damping  ratio  and  03   represents  the 

d 
desired  natural  frequency  of  the  model  system.   The  output 

of  the  model  is 


cm(t)  =  Yl(t)  .  (III. 4) 

A  performance  index  which  is  suitable  for  many  design 
problems  is 


J(a)  = 


r   f  2 

(c  (t).  -  c  (t))  dt  ,  (III. 5) 

O  III 

o 


where  ct.   is  an  N  -dimensional  vector  representing  the  free 
system  parameters,  and  c  (t)  represents  the  model  response 
(for  the  case  under  discussion  it  is  the  output  of  the 
second-order  model).  .  t_  (the  final  time  of  the  integra- 
tion) can  be  taken  as  the  settling  time  of  the  second- 
order  model;   it  can  be  calculated  by  using  the  equation 

fcf  =  ri-  =  ^  (III-6) 

d  n .,    d 


where  cr,  is  the  real  part  of  the  desired  root  locations. 

If  another  model  response  is  used,  the  final  time  can 
be  guessed  by  using  the  performance  specifications  —  it 
should  be  long  enough  to  cover  at  least  one  period  of  the 
lowest  frequency  of  the  transients  which  occur  during  the 
minimization  process.   The  value  of  the  performance  index 
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depends  on  tf,  but  once  it  is  chosen  (long  enough)  and 
fixed,  it  does  not  affect  the  optimum  parameter  values  ob- 
tained at  the  end  of  the  minimization.   If  tf  is  chosen 
longer  than* is  necessary,  extra  computer  time  may  be  spent. 
For  this  reason,  if  the  designer  does  not  have  any  idea 
about  the  dynamic  behavior  of  the  system,  he  should  simulate 
the  system  by  using  the  starting  values  of  the  free  system 
parameters  as  the  first  step  of  the  design  procedure. 

The  performance  index  in  the  sampled  form  is  written 
as 

nf 

J (Sf)  =  Y   (cs(kT>  "  cm(kT))2  (ill. 7) 

k=0 


where  nf  is  a  positive  integer  defined  a 


s 


final  time  of  the  calculations     f       i-r-r-r    o\ 

nc  =   5 =-: ; — -5 =  —  .     (III. 8) 

f         Sampling  period  T 


When  J  is  minimized  with  respect  to  the  free  system 
parameters  (a)  in  the  bounded  region  defined  by  the  con- 
straints (see  Section  III.C),  the  optimum  parameter  values 
(a*)  are  found. 

The  form  of  the  performance  index  can  be  modified  to 
fit  the  needs  of  the  specific  design  problem;   variation 
from  the  basic  form  (Eqs.  (III. 5)  and  (III. 7)  are  employed 
in  the  examples  given  in  the  later  parts  of  the  thesis 
(see  particularly  Examples  7  and  9)  . 
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E.    EXAMPLE  5 

The  linear  system  design  problem  solved  as  Example  1 
in  Chapter  II  is  reconsidered  here.   This  time,  a  solution 
is  obtained  by  using  the  method  "Optimization  for  the  Best 
Response  in  the  Time-Domain".   The  system  (shown  in  Fig. 
11.3(a))  and  all  specifications  are  kept  the  same  to  pro- 
vide a  basis  for  comparison.   A  step-by-step  solution  is 
given  below; 

1.   The  Model  Response 

The  desired  dynamic  response  was  represented  by 
two  complex  roots,  in  Example  1;   the  desired  dynamic  re- 
sponse (the  model  response)  is  the  output  of  a  second-order 
model  system  defined  by  dominant  roots  located  at  cr  =  0.5, 
U),  =  ±  0.5.   The  state  equations  of  the  second-order  model 
are  written  from  Eq.  (III. 3) 


x  (t)  =  x  (t) 
b  b  (III. 10) 

x6(t)  =  0.5(r(t)  -  x5(t))  -  x6(t)  , 

where  x,-(t)  =  y-i(t)  and  x,  (t)  =  y~(t)  are  used  for  program- 
ming convenience.   A  unit  step  input  is  assumed  (r(t)=  l(t)) 
The  model  response  is 


c  (t)  =  xc(t)  .  (III. 11) 

m       o 


2.   Type  of  Compensation 

The  same  cascade  compensator  (shown  in  Fig. II. 4 (a)) 
is  used  for  compensation. 
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3 .   State  Equations  of  the  System 

The  state  equations  of  the  system  are  written  from 
the  transmission  function  (Eq.  (11.11)  in  the  rational  ca- 
nonical form 

xx(t)  =  x2(t) 

x9(t)  =  x  (t) 

2       J  (III. 12) 

x3(t)  =  x4(t)  +  h3  r(t) 

x4(t)  =  -bQx1(t)-  blX2(t)-  b2x3(t)-  b3x4(t)+  h4r(t) 
where 


bQ  =  Kp  b3  =  2  +  p 

h1   =   p  +  K/k  h3  =  K/k  (III. 13) 

b2  =  1  +  2p  h4  =  bQ  -  h3b3  . 


The  output  of  the  system  is 

c  (t)  =  x,  (t)  .  (III. 14) 

s       l 

4.   The  Performance  Index 

The  performance  index  defined  in  Eq.  (III. 5)  is 
suitable  for  this  design  problem,  hence 

fcf 
J  £  x7(t)  =  J   (xx(t)  -  x5(t))2dt  ,  (III. 15) 

and  its  time  derivative  is 

x?(t)  =  (xx(t)  -  x5(t))2  .  (III. 16) 
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The  final  time  for  the  integration  is  taken  as  the 
settling  time  of  the  model  response 

tf  =  4/cr   =  8  seconds.  (III.  17) 

5 .  Constraints 

The  only  constraining  inequalities  are 

0.1  ^  k  <  10 

p  >  0  (III. 18) 

0.4  ^  K  <  0.5  . 

6.  Minimization 

The  performance  index  (x_(t))  is  minimized  with 
respect  to  the  free  system  parameters  (a) .   A  gradient 
method  was  used  for  the  minimization  (COMPUTER  PROGRAM  II, 
given  at  the  end  of  the  thesis  is  typical  for  this  method) . 
At  each  step  of  the  minimization  procedure,  Eqs.  (III. 10), 

(III. 12)  and  (III. 16)  are  integrated  from  t  =  0  to  'tf  by 

1 
using  a  fourth-order  Runge-Kutta  method  .   The  optimum 

parameter  values  are  k*  =  0.1,  p*  =  9.84,  and  K*  =  0.4945. 

The  optimum  time  response  of  the  system  and  the  desired 

response  are  shown  in  Fig.  III.l. 

The  system's  response  shown  in  this  figure  (and  its 

closeness  to  the  desired  response)  should  be  compared  with 

the  response  shown  in  Fig.  II. 7. 


SUBROUTINE  RKLDEQ,  N.P.G.S.  Computer  Facitlity  Program 
Library. 


56 


L.    t . 

Eh 

p 

En 
P 

0.8        ° 

— — — 

1 

/ 

_.     j 

<*K 

0.6 

/ 

0.A _ / 

-4 

I 

1 

/' 

r 

i 

f 

! 
1 

1 

0,2    J_ 
I 

I  J 

2. 

-  — 

4. 

. 

6. 

10. 

8. 

Fig.  III.l.  EXAMPLE  5 


Compensated  system's  response  (A) 

and  the  desired  response  (B)  to  a 
unit  step  input. 
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F.    EXAMPLE  6 

Example  2  of  Chapter  II  is  re-solved  by  using  the 
method  "Optimization  for  the  Best  Response  in  the  Time- 
Domain".   The  system  is  shown  in  Fig.  II. 8. 

In  the  solution  given  below,  paragraph  numbers  indic- 
ate the  corresponding  steps  in  Example  5. 

1.   The  model  response  is  obtained  from  the  output  of 
a  second-order  model  system  which  is  defined  by  the  state 
equations 


x5(t)  =  x6(t) 

x6(t)  =  5(r(t)  -  x5(t)  -  2x6(t)  , 


(III. 19) 


and  the  model  response  is 

c  (t)  =  x_(t)  .  (III. 20) 

m       o 

A  unit  step  input  is  used  in  the  design,  thus 

r(t)  =  i(t)     .  (III. 21) 

2.  The  same  tachometer  and  acceleration  feedbacks  are 
used  to  compensate  the  system. 

3.  The  state  equations  of  the  system  are  written 
directly  from  Fig.  11.8(b). 

x^t)  =  x2(t) 

X  (t)   =  X   (t) 

2  J  (III. 22) 

x3(t)  =  x4(t) 

x4(t)  =  [r(t)-x1(t)-Ktx2(t)-Kax3(t)]K-[50x2(t)+65x3(t)  +  16x4(t)] 
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The  output  of  the  system  is 


cs(t)  =  x1(t)  .  (III. 23) 

4.   The  time  derivative  of  the  performance  index  is 

x7(t)  =  (x1(t)  -  x5(t))2  .  (III. 24) 
The  final  time  for  the  integration  is 


t   =  4/cr   =  4/1  =  4  seconds.  (III. 25) 

5.   The  constraints  on  the  free  variables  are 


0  ^  K   ^  10  , 


0  £  K  <■   10  , 
a 


(III. 26) 


and  the  system  gain  is  set  to  its  desired  value  (K  =  1000) . 
6.   The  performance  index  defined  by  the  integral  of 

Eq.  (III. 24)  is  minimized  by  using  a  gradient  method;   the 

* 

optimum  parameter  values  are  K*  =  0.142,  K  =  0.4,  and  the 

a  t 

optimum  time  response  of  the  system  (for  a  unit  step  input) 
is  plotted  in  Fig.  III. 2.   In  the  same  figure  the  desired 
response  is  also  shown.   The  result  of  s-domain  synthesis 
procedure  was  given  in  Fig„  II. 9. 

G.    EXAMPLE  7 

The  linear  system  design  problem  solved  in  Examples 
2,  3,  and  6  is  re-considered  here.   As  in  Example  3  of 
Chapter  II  dynamic  performance  specifications  require  "the 
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Fig.  III. 2.  EXAMPLE  6.   Compensated  system's  response  (A)  and 

the  desired  response  (B)  to  a  unit  step 
input. 
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fastest  response"  to  a  unit  step  input.   This  problem  is 
solved  three  times  by  using  slightly  different  techniques 
to  demonstrate  the  effects  of  the  different  forms  of  per- 
formance indices. 
Solution  1 

1.  .The  unit  step  input  is  used  as  the  model  response: 

cm(t)  =  r(t)  =  £(.t)  .  (III. 27) 

Other  parts  of  the  solution  are  exactly  the  same  as 
the  one  given  in  the  previous  example  except  that  the  per- 
formance index  now  takes  the  form 

x7(t)  =  (l.-x1(t))2  (III. 28) 

2.  The  new  performance  index  is  minimized  with  respect 

to  the  free  system  parameters  (K   and  K  )  ant  the  optimum 

a       t 

parameter  values  are  K*  =  0.0296,  K*  =  0.4546.   The  optimum 

a  t 

response  is  plotted  in  Fig.  III. 3  .   In  the  same  figure  the 
desired  response  of  Example  6  is  also  shown,  for  comparison. 

Solution  2 

The  optimum  parameter  values  in  Solution  1  yield  a 
faster  response  (shorter  rise  time  and  shorter  settling 
time)  than  obtained  in  Example  6,  but  the  response  is  some- 
what oscillatory  and  the  settling  time  is  about  four  sec- 
conds.   If  dynamic  performance  specifications  require  a  less 
oscillatory  output  and  a  settling  time  less  than,  for  exam- 
ple, 2.5  seconds,  a  small  change  in  the  design  procedure 
yields  the  desired  dynamic  response. 
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Fig.  I II. 3.  EXAMPLE  7 


The  fastest  time  response  (A)  and  the 
desired  response  (B)  (for  EXAMPLE  6) 
to  a  unit  step  input.  (SOLUTION  1). 
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For  this  application  the  performance  index  takes  the 
form 


J  =  x?(t)  = 


r   f  2 

w  (l.-x1(t))  dt  ,  (III. 29) 

o 


where  tf  =  4  seconds  as  before  and  w.  is  the  penalty  func- 
tion defined  as 

for   t  <  2.5 
w1      =  1  (III. 30) 

1000   for   t  ^  2.5  . 

This  penalty  function  assures  a  settling  time  less  than  or 

equal  to  2.5  seconds.   Minimization  of  Eq.  (III. 29)  yields 

the  optimum  parameter  values;   K*  =  0.047,  K*  =  0.441. 

at 

The  output  of  the  compensated  system  with  optimum  parameter 
values  is  shown  in  Fig.  III. 4. 

Solution  3 

A  somewhat  different  approach  to  the  same  design  prob- 
lem is  considered.   The  differences  from  the  previous 
solutions  are  given  below. 

The  system  gain  is  used  as  the  third  free  parameter 

3         4 
(which  is  constrained  in  the  region  10   ^  K  ^  10  )  . 

The  performance  index 

J  =  x?(t)  =  |   [w1(l.-x1(t))2  -  x^(t)  |  dt      (III. 31) 
1  o 

is  used  for  the  design.   In  this  approach  the  time  deriva- 
tive of  the  output  is  maximized  to  obtain  the  fastest 
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Fig.  III. 4.  EXAMPLE  7. 


The  fastest  time  response  (A)  (with  a 
constraint  on  the  settling  time)  and 
the  desired  response  (B)  (for  EXAM- 
PLE 6)  to  a  unit  step  input. 
(SOLUTION  2) . 
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2 

response.   The  first  term  (w  (l.-x,(t))  )  is  used  to  satis- 
fy another  (assumed)  performance  specification  which  requires 
the  peak  overshoot  be  less  than  or  equal  to  40%;   for  this 
purpose  the  penalty  function 


w 


for   I  xx(t)  I  ^  1.4 

1     )  for  all  t  in  [0, t  ] 

100    for   |  x1(t)  |  >  1.4 

(III. 32) 

is  used. 

Minimization  of  the  performance  index  (Eq.  (III. 31)) 

yields  the  optimum  parameter  values  K*  =  1012.6,  K*=0.07, 

a 

K?  =  0.296.   The  optimum  response  is  shown  in  Fig.  III. 5. 

H.    EXAMPLE  8 

In  the  Example  4  of  Chapter  II  a  third-order  linear 
system  was  compensated  for  the  given  dominant  root  loca- 
tions.  The  use  of  a  lag  compensator  (in  the  s-domain)  was 
demonstrated  in  Example  4.   The  same  example  is  re-solved 
here  to  compare  the  results  of  the  two  methods. 

The  system  is  shown  in  Fig.  II. 11(a).   Paragraph 
numbers  are  used  to  indicate  the  corresponding  steps  in 
Example  5 . 

1.   The  desired  root  locations  in  Example  4  (cy      =   1, 
to,  =  ±  2)  correspond  to  a  second-order  system  which  is 
defined  by  the  state  equations 

x  (t)  =  x  (t) 

b       b  (III. 32) 

x6(t)  =  5(r(t)  -  x5(t))  -  2x6(t)  . 
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Fig.  III. 5.  EXAMPLE  7 


The  fastest  response  (A)  by  maximizing 
the  derivative  term,  and  the  desired 
response  (B)  (for  EXAMPLE  6) 
(SOLUTION  3) . 
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The  model  response  is 


c  (t)  =  x  (t)  . 
m        b  . 


(III. 33) 


The  unit  step  input  is  used  in  the  design,  i.e. 


r(t)  =  i(t)     . 


(III. 34) 


2.  -  3.   A  cascade  compensator  is  used  (as  shown  in 
Fig.  11.11(b))  to  obtain  the  desired  response.   The  state 
equations  of  the  system  are  written  from  the  system's 
transmission  function  (Eq.  (11.45))  in  the  rational  canon- 
ical form 


x(t)  =  Ax(t)  +  br(t) 


(III. 35) 


where 


A   = 


0  10  0 
0  0  10 
0    0    0    1 


•cl    — cln    a. «    a.  ~ 

o    1    2    3 


(III. 36) 


b  =  (0  0  h_  h„)T  ,  and 


a   =  Kp 


a3=  12  +  p 


a   =  36p+K/k   h3=  K/k 


a^  = 


36+12p   h„=  a  -  h0a^ 
^  4    o    3  3 


(III. 37) 


The  output  of  the  system  is 


cg(t)  =  xx(t) 
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4.  The  time  derivative  of  the  performance  index  is 
x7(t)  =  (x1(t)  -  x5(t))2  .  (III. 38) 

5.  The  constraints  on  the  free  system  variables  are 

425  £  K  £  1000 

0.1  £  k  ^  10  (III. 39) 

0<p  f   100  (arbitrary)  . 

6.  The  performance  index  defined  by  the  integral  of 
Eq. (III. 38)  is  minimized  with  respect  to  the  free  system 
parameters  and  the  optimum  parameter  values  are  K*  =  554, 
k*  =  7.84,  and  p*  =  0.0116.   The  optimum  response  to  a 
unit  step  input  and  the  model  response  are  shown  in  Fig. 
III. 6. 

I.    CONCLUDING  REMARKS 

In  this  chapter,  the  method  "Optimization  for  the  Best 
Response  in  the  Time-Domain"  has  been  developed  and  applied 
to  several  linear   system  design  problems.   The  same  exam- 
ples, given  in  Chapter  II  for  the  s-domain  design  method, 
have  been  reconsidered  here  to  make  comparison  of  the  two 
design  methods  possible.   From  this  comparison  and  from 
the  general  development  of  the  method,  it  may  be  concluded 
that  this  method  has  many  advantages  over  the  present  dynam- 
ic system  design  procedures.   The  most  important  advantages 
are  listed  below: 
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Fig.  III. 6.  EXAMPLE  8. 


Time  response  (A)  of  the  compensated 
system  with  a  lag  compensator  and  the 
desired  response  (B) ,  to  a  unit  step 
input. 
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1.  This  method  is  applicable  to  high-order,  multi- 
parameter systems.   Any  dynamic  system  which  can  be  repre- 
sented by  state  equations  can  be  designed  by  this  method; 
the  order  of  the  system  and  the  number  of  free  system  para- 
meters can  be  very  high  without  introducing  any  difficulty 
in  the  design  procedure. 

2.  The  work  required  to  design  a  high-order  system  is 
less  than  with  other  methods;   in  fact,  if  a  computer  pro- 
gram for  the  minimization  process  is  available,  the  only 
required  work  is  in  obtaining  the  state  equations  of  the 
system. 

3.  This  method  directly  shapes  the  time  response  of 
the  system  by  adjusting  the  free  parameters.   If  the  op- 
timum response  obtained  by  this  method  is  not  close  enough 
to  the  desired  (model)  response,  then  the  system  configura- 
tion must  be  changed.   The  only  indirect  part  of  the  method 
is  that  the  designer  must  decide  on  the  available  system 
components  and  compensation  techniques. 

The  method  "Optimization  for  the  Best  Response  in  the 
Time-Domain"  enjoys  the  following  advantages  over  the  s- 
domain  method  developed  in  Chapter  II. 

a.  The  basic  idea  behind  this  method  is  much  sim- 
pler. 

b.  Optimization  for  the  best  response  in  the  time 
domain  yields  a  better  answer  to  the  design  problem  because 
it  is  more  direct  and  it  does  not  employ  any  questionable 
assumptions  such  as  dominancy  of  the  desired  roots. 
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c.  This  method  requires  less  hand  calculation  and 
less  programming  effort. 

d.  This  method  is  not  limited  to  the  design  of 
linear  systems.   Besides  its  simplicity  this  is  the  most 
important  advantage  of  the  method,  because  it  provides  a 
powerful  computerized  approach  to  the  design  of  nonlinear 
and  sampled-data  systems  and  also,  for  example,  systems  with 
time  delay.   These  applications  are  presented  in  Chapter  IV. 
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IV.   OPTIMIZATION  FOR  THE  BEST  RESPONSE  IN  THE 
TIME  DOMAIN  —  APPLICATIONS  TO  NONLINEAR 
AND  SAMPLED-DATA  SYSTEMS  AND  SYSTEMS 

WITH  TIME  DELAY 

, r, 

A .  GENERAL 

The  method  "Optimization  for  the  Best  Response  in  the 
Time-Domain"  is  applicable  to  a  large  variety  of"  dynamic 
system  synthesis  problems.   Any  dynamic  system  which  can 
be  represented  by  state  equations  can  be  designed  to  yield 
a  desirable  dynamic  and  steady-state  performance. 

In  this  chapter,  applications  of  the  method  to  non- 
linear and  sampled  data  systems  and  to  systems  with  time 
delay  are  considered.   Three  simple,  but  representative, 
examples  are  given.   Application  to  these  classes  of  sys- 
tems increases  the  value  of  the  method,  because  direct  and 
efficient  design  methods  are  rare  in  these  areas  of  control 
theory. 

If  the  system  is  not  linear,  it  may  require  some  addi- 
tional logic  (in  the  computer  programming)  to  obtain  the 
state  values  in  the  time  interval  of  interest;   but  the 
main  procedure  (as  presented  in  Chapter  III)  is  directly 
applicable. 

B.  EXAMPLE  9 

A  chemical  process  with  an  inherent  time  delay  (T  =  0.5 
seconds)  is  to  be  regulated.   The  system  is  shown  in  Fig, 
IV. 1.   The  load  changes  are  sudden  and  can  be  assumed  to  be 
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step  functions.   The  only  available  free  parameter  is  the 
amplifier  gain  (K  ) .   The  desired  dynamic  and  steady-state 
specifications  will  be  defined  later. 

For  the  given  system  configuration  the  system's  re- 
sponse can  be  divided  into  two  parts: 

1.  Steady-State  Accuracy.   Since  the  given  system  is 
type-O,  there  will  always  be  a  constant  steady-state  error 
for  step  disturbances.   The  value  of  the  steady-state 
error  (E   ) ,  which  depends  on  the  amplifier  gain  (K  ) ,  is 

SS  C 

E   (K  )  =  1/(1+4.8  K  )  (IV. 1) 

ss   c     '         c 

2.  Dynamic  Accuracy.   This  is  also  a  function  of  the 
amplifier  gain  (K  ) .   The  dependence  of  the  shape  of  the 
dynamic  response  (on  K  )  will  be  apparent  at  the  end  of  the 
solution. 

In  general,  high  gain  values  cause  more  oscillations, 
and  less  steady-state  error  (see  Eq.  (IV.l));   low  gain 
values  result  in  a  less  oscillatory  response  but  more 
steady-state  error.   With  only  one  free  parameter,  it  is 
impossible  to  meet  all  of  the  specifications;   therefore, 
the  design  method  should  seek  the  best  compromise  between 
the  desired  specifications. 

To  make  the  initial  conditions  of  the  states  all  zero, 
it  can  be  assumed  that  the  reference  signal  is  zero  (u(t)=0)), 
and  that  the  system  is  in  a  steady-state  condition.   The 
system's  behavior  in  case  of  sudden  load  changes  will  be  in- 
vestigated. 
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Fig.  IV. 1.  EXAMPLE  9 


The  state  equations  of  the  system  are 


x^t)  =  x2(t)+  r(t-0.5)/3 


(IV. 2) 


x2(t)  =  -(9.6Kc  x1(t-0.5)+2x1(t)+9x2(t)+r(t-0.5))/9 

There  are  three  terms  affected  by  the  0.5  -  second 
time  delay.   x, (t-0.5)  represents  the  value  of  x,(t),  0.5 
seconds  earlier  during  the  process.   The  load  disturbances 
are  assumed  to  be  unit  step  functions;   therefore, 


r(t-0.5)  = 


0   for   t  <  0.5 


1   for   t  *  0.5 


(IV. 3) 


can  be  used. 

The  output  of  the  system  is 


cs(t)  =  x1(t)  . 


(IV. 4) 


It  the  performance  specifications  only  require  that 
the  output  of  the  system  should  be  kept  as  close  to  the 
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reference  as  possible,  and  the  transients  should  last  not 
more  than  30  seconds,  then  a  performance  index,  which  meas- 
ures the  deviation  of  the  output  from  the  reference  level 
(from  t  =  0  to  t.p=  30  seconds)  can  be  defined  as 

c'     3 


^x^(t)  ^     (x1(t)  -  u(t)rdt  .         (IV. 5) 
o 


Since  the  reference  level  was  assumed  to  be  zero,  this  per- 
formance index  becomes 


x3(t)  = 


?    f   2 

x],  (t)  dt  .  (IV. 6) 

o 


The  time  derivative  of  the  performance  index  is 

x3(t)  =  X;L2(t)  .  (IV. 7) 


Minimization  of  this  performance  index  with  respect  to 
K  yields  the  "best"  response  of  the  system  within  the 
limitations  of  the  system  configuration. 

Eqs.  (IV. 2)  and  (IV. 7)  are  integrated  with  respect  to 
time  to  calculate  the  value  of  the  performance  index  at 
each  step  of  minimization  process.   The  integration  step 
size  is  chosen  as  an  integer  submultiple  of  the  time  in  the 
process  (At  =  0.05  is  chosen).   For  the  first  ten  integra- 
tion steps  the  values  of  r (k  At  -  0.5)  and  x, (k  At  -  0.5) 
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are  zero  ;   the  value  of  the  variable  x, (k  At)  are  stored 
(k  is  the  integrating  step  counter) .   After  the  tenth  in- 
tegration step  (k  >  10),  r (k  At  -  0.5)  is  equated  to  unity 
and  the  values  of  x,  ((k-10)At)  are  found  and  used  from  the 
stored  values  of  x, (k  At) . 

The  only   difference  from  a  standard  linear  system  de- 
sign problem  is  in  the  logic  of  expressing  the  existing 
time  delay  in  the  system,  and  as  explained  above,  this  does 
not  introduce  any  complication  to  the  solution. 

K*  =  1.66  is  found  by  the  minimization  procedure.   The 
c 

time  response  of  the  system  to  a  unit  step  load  change  is 
shown  in  Fig.  IV. 2.   Peak  overshoot  is  about  33%  and  the 
steady-state  error  is  11.1%. 

If  the  desired  response  is  to  be  less  oscillatory,  one 
can  include  the  time  derivative  of  the  output  in  the  perform- 
ance index;   that  is, 

x3(t)  =  X;L2(t)  +  xx2(t)  ,  (IV. 8) 

or  by  substituting  the  value  of  x,  (t) (from  Eq.  (IV. 2))  one 
obtains 

x3(t)  =  xx2(t)  +  (x22(t)  +  r(t  -  0.5)/3)2  .    (IV. 9) 


Since  the  numerical  integration  is  a  discrete  process, 
the  arguments  of  the  delayed  variables,  r(t-0.5)  and 
x, (t-0.5)  are  represented  in  discrete  form. 
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Fig.  IV. 2.  EXAMPLE  9.  Time  response  of  the  system  to 

a  unit  step  disturbance  (j  =  xj) 
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By  minimizing  this  performance  index  K*  =  1.5  is  ob- 

c 

tained  (The  time  response  is  shown  in  Fig.  IV. 3).   This 
time  response  is  less  oscillatory  but  suffers  more  steady- 
state  error  (12.2%). 

To  meet  the  given  specifications,  the  terms  represent- 
ing the  output  and  its  time  derivative  (in  the  performance 
index)  may  be  multiplied  by  suitable  weighting  factors. 
For  example,  K*  =  1.12  and  the  time  response  shown  in  Fig. 
IV. 4  are  obtained  by  minimizing  the  performance  index 

x3(t)  =  w1  xx2(t)  +  w2  xx2(t)  ,  (IV. 10) 

for  w,  =  1   and  w„  =  10. 

If  the  desired  response  should  have  the  smallest  cu- 
mulative rate  of  change  (slowest  varying  transients  in  the 
time  interval  [o,tf]),  the  performance  index  takes  the  form 

x3(t)  =  xx(t)  =  (x2(t)  +  r(t  -  0.5)/3)2  .      (IV. 11) 

This  performance  index  is  the  same  as  (IV. 10),  of  course, 
when  w   =  0  and  w  =  1  is  used.   With  this  performance  in- 
dex K*  =  0.473  is  obtained  (for  the  time  response  see  Fig. 
IV. 5). 

Finally  some  applications  may  require  that  the  system's 
output  should  not  deviate  from  the  reference  level  more  than 
a  certain  percentage  for  a  given  load  change.   To  satisfy 
this  requirement  a  constraint  on  the  output  is  inserted. 
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Fig.  IV. 3.  EXAMPLE  9.  Time  response  of  the  system  to  a  unit 

step  disturbance  (J  =  x.   +  x,  ) . 
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Fig.  IV. 4.  EXAMPLE  9.  Time  response  of  the  system  to  a 

unit  step  disturbance  (J  =  x,+10x,) 
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Fig.  IV. 5.  EXAMPLE  9.  Time  response  of  the  system  to  a 
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For  example,  to  ensure  that  the  maximum  deviation  from  the 
reference  level  is  less  than  31%  of  the  load  disturbance, 
the  performance  index  defined  in  Eq.  (IV. 8)  is  used  and  the 
additional  constraint 

|x1(t)  |  ^  0.31   for  all   t   in  [0,tf]         (IV. 12) 

.  .  .    .   1   .       * 
is  imposed.   Minimization  yields  K  =  1.85;   the  time  re- 
sponse of  the  system  (for  a  unit  step  load  change) is  shown 
in  Fig.  IV. 6. 

This  series  of  solutions  (see  Table  IV. 1)  shows  that 
a  system  can  be  optimized  for  any  desired  response  (that  is, 
the  most  suitable  response  within  the  limitations  of  the 
given  system  configuration  can  be  obtained)  by  using  the 
method  "Optimization  for  the  Best  Response  in  the  Time- 
Domain"  —  the  existence  of  one  or  more  time  delays  does 
not  introduce  any  difficulty  or  complication  in  the  solu- 
tion process. 


The  penalty  function  technique  used  in  Example  7  was 
applied. 
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J 

K* 

c 

Peak 
overshoot 

E 
ss 

Time 
response 

2 
Xl 

1.66 

33% 

11.1% 

Fig.  IV. 2 

2    •  2 
xl   +  Xl 

1.5 

35% 

12.2% 

Fig.  IV. 3 

X]L2  +  10  xx2 

1.12 

36% 

15.6% 

Fig.  IV. 4 

.  2 
Xl 

0.473 

45% 

30% 

Fig.  IV. 5 

2    .  2 
Xl   +  Xl 

1.85 

31% 

10  .  1% 

Fig.  IV. 6 

1  x   £  31% 

for  all  t 

in  the  time 

interval 

[o,tf] 

TABLE  IV. 1.  EXAMPLE  9  (SOLUTIONS) 

C.    EXAMPLE  10 

A  simple  second-order  sampled-data  system  (shown  in 
Fig.  IV. 7 (a))  is  considered.   The  two  adjustable  design 
parameters  (free  system  parameters) are  the  system  gain  (K) 
and  the  sampling  period   (T) . 


1   The  method  "Optimization  for  the  Best  Response  in  the 
Time-Domain"  yields  an  important  advantage  when  it  is  ap- 
plied to  the  sampledrdata  systems  —  the  sampling  period 
can  be  taken  as  one  of  the  free  variables  without  increas- 
ing the  complexity  of  the  solution. 
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Fig.  IV. 7 (a).  EXAMPLE  10.   The  system, 


Fig.  IV.7(b). 

EXAMPLE  10. 

Time  response  of  the 
system  (A)  without 
a  hold  circuit  and 
the  desired  response 
(B).  , 
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Solution  (without  a  Hold  Circuit) 

2 

1.   The  dynamic  performance  specifications   require 

that  the  system's  response  should  be  as  close  to  the  re- 
sponse of  the  second-order  continuous  system  as  possible. 
The  state  equations  of  the  second-order   continuous  system 
are  ' 


x^(t)  =  x  (t) 

(IV. 13) 

x4(t)  =  0.526(r(t)  -  x3(t))  -  0.630x4(t)  , 
and  the  model  response  is 


c  (t)  =  x_(t)  .  (IV. 14) 

m       J 

The  unit  step  input  is  used  in  the  design,  i.e. 

r(t)  =  i(t)  .  (IV. 15) 

3.   The  state  equations  of  the  system  are  written  from 
Fig.  IV. 7(a) 

xx(t)  =  x2(t) 

(IV. 16) 
x2(t)  =  K  e^t)  -  x2(t)  . 

Since  the  hold  circuit  is  not  used  in  this  solution 
(see  Fig.  IV. 7 (a) ) . 


Paragraph  numbers  indicate  the  corresponding  steps  in 
Example  5. 
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e  (t)  =  e2(t)  ,  (IV. 17) 

and  also 

fO       if  the  sampling  switch  is  open 
ex(t)  =  ) 

(r(t)-x,(t)   during  the  sampling  instant 

(IV. 18) 
Logic  must  be  included,  in  the  computer  program,  to 
find  the  sampling  instants,  during  the  integration;   in 
fact  this  is  the  only  required  addition  to  the  computer 
program  which  is  written  for  a  standard  linear  system  de- 
sign. 

It  is  assumed  that  the  sampling  switch  is  closed  for 
the  first  time  at  t  =  0,  and  opened  after  an  integration 
step.   In  other  words  at  the  first  step  of  the  discrete 
numerical  integration  process  e,  (t)  =  r(t)  -  x.  (,fe.)  is  used, 
then  r, (t)  is  set  to  zero.   Numerical  integration  is  con- 
tinued for  "I  -  1"  steps  with  the  e  (t)  =0  value.   I  is  an 
integer  found  from  the  equation 

I  =  Sampling  Period /A  (IV.19) 

Integration  step  size 

A  step  counter  (i)  counts  the  integration  steps  starting 
from  the  sampling  instant.  i=  I  -  1  indicates  another  sam- 
pling instant v  when  this  occurs 


i  =  0  ,  and 

(IV. 20) 
ex(t)  =  e(t)  -  x1(t) 
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are  set  by  the  computer  program,  and  a  new  count  is  started. 
The  sampling  time  is  assumed  to  be  one  integration  step  size 
in  this  solution.   If  the  sampling  time  (the  period  of  time 
that  the  sampling  switch  stays  closed)  should  be  long  com- 
pared to  the  integration  step  size  then  another  counter  may 
be  added  to  indicate  the  opening  instants  of  the  sampling 
switch.   The  integration  step  size  should  be  chosen  as  an 
integer  submultiple  of  both  the  sampling  period  and  sampling 
time. 

The  output  of  the  system  is 


cg(t)  =  xx(t)  .  (IV. 21) 


4.  The  time  derivative  of  the  performance  index  is 

x5(t)  =  (x-^t)  -  x3(t))2  .  (IV. 22) 

The  final  time  of  the  integration  is  taken  as  the  set- 
tling time  of  the  continuous  model, 

t   =  4/(7   =  12.7  seconds  .  (IV. 23) 

5.  Constraints  on  the  free  variables  are 

0.1  £  T  ^  1.5 

(IV. 24) 
0  <  K  £  15  . 

6.  The  performance  index  defined  by  the  integral  of 
Eq.  (IV. 22)  is  minimized  with  respect  to  the  free  system 
parameters  (T  and  K)  and  the  optimum  parameter  values  are 
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T*  =  0.6  ,   K*  =  7.20.    The  output  of  the  system  to  a 
unit  step  input  and  the  model  response  are  shown  in  Fig. 
IV. 7(b) . 

Solution  (With  a  Hold  Circuit) 

If  a  zero-order  hold  circuit  is  inserted  after  the 
sampling  switch  the  only  necessary  change  is  in  the  logic 
given  in  conjunction  with  the  Eq.  (IV. 18). 

After  a  sampling  instant,  e~(t)  is  not  equated  to 
zero,  but  it  is  retained  at  the  value  of  the  last  sampling 
for  the  remaining  I  -  1  steps  of  integration.   In  this 
solution  impulse  sampling  is  used,  i.e.,  the  switch  is  not 
kept  closed  for  0.05  seconds  but  it  is  opened  as  soon  as 
the  sample  is  measured. 

Minimization  of  the  same  performance  index  (Eq. (IV. 22)) 
yields  T*  =  1.05  and  K*  =  0.5875.   The  output  of  the  sam- 
pled data  system  with  zero-order  hold  and  the  model  response 
are  shown  in  Fig.  IV. 8. 

D.    EXAMPLE  11 

A  designer  simulates  the  fourth-order  linear  system 
shown  in  Fig.  IV. 9  and  obtains  the  time  response  to  a  unit 
step  input  shown  in  Fig.  IV. 10.   This  response  is  for  a 
system  gain  of  K  =  1600.   It  is  assumed  that  this  dynamic 
response  perfectly  fits  the  particular  application, except 
that  the  velocity  coefficient  must  be  doubled  for  steady- 
state  accuracy.  When  the  gain  is  doubled  (K  =  3200)  however, 
the  system  moves  to  its  stability  limit,  and  the  time  re- 
sponse is  quite  oscillatory  as  shown  in  Fig.  IV. 10. 
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Fig.  IV. 8.  EXAMPLE  10. 


Time  response  of  the 
system  (A)  with  a  zero 
order  hold  circuit, 
and  the  desired 
response  (B) . 
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Fig.  IV. 9.  EXAMPLE  11.  The  linear  system. 

The  real  system  has  two  inherent  nonlinearities  —  a 
saturating  amplifier  in  the  forward  path  and  dead  zone  in 
the  feedback  path.   The  real  system  with  the  nonlinearities 
is  shown  in  Fig.  IV. 11.   The  designer  decides  to  use  two 
stages  of  linear  cascade  compensators  to  obtain  a  time  re- 
sponse as  close  to  the  response  of  the  linear  system  as 
possible;   in  other  words,  the  linear  system  shown  in  Fig. 
IV. 9  (with  K  =  1600)  will  be  used  to  generate  "the  model 
response"  for  the  optimization  process. 

The  compensated  system  has  four  free  parameters  (see 
Fig.  IV. 12)  —  k, ,  k~  (the  compensation  ratios  of  the  first 
and  second  compensators)  and  p, ,  p~  (the  pole  locations  of 
the  compensators) . 

The  state  equations  of  the  ideal  (linear)  system  are 


x?(t) 


x8(t) 


x8(t) 


*9(t) 


x9(t) 


x1Q(t)  +  K±    (r(t)  -  x?(t)) 


(IV. 25) 


x.0(t)  =  -(lOOOx  (t)+350x  (t)+35x1  (t)+27Ki(r(t)-x7(t)) 
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Fig. IV. 10.  EXAMPLE  11 


Time  response  of  the  nonlinear  system 
(A)  (for  K  =  3200),  and  the  time  re- 
sponse of  the  linear  system  (B)  (for 
K  =  1600) ,  (unit  step  input  is  used) . 
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r(t)   _y2(t) 


+ 


z2(t) 


x1(t) 


xx(t) 


where  the  PLANT  is  defined  by  the  transfert  function 


Z1(s) 


=  G(s)  = 


s+8 


s(s+5) (s+10) (s+20) 


(a)   Real  (Nonlinear)  system 


►  Vo(t) 


z2(t) 


-K 


K  =  3200,  K  =  1600,  K  =-1400      K   =  0.2 
x  n  c 


slope  =  1 


K 


'  x]_(t) 


(b)   N   -  Block 


(c)   N2  -  Block 


Fig.  IV. 11.   EXAMPLE  11.   The  real  system. 
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where  K.  =  1600  and  r(t)  =  ll(t)  —  the  unit  step  function. 
The  model  response  is  the  output  of  the  linear  system 


c  (t)  =  x  (t)  .  (IV. 26) 

m       / 

The  state  equations  of  the  real  (nonlinear)  system 
with  two  cascade  compensators  (refer  to  Fig.  IV. 12)  are 
written  in  the  following  steps: 

1.   The  plant  is  specified  by 


x1(t)  =  x2(t) 


x2(t)  =  x3(t) 

x3(t)  =  x4(t)  +  Z;L(t) 


(IV. 27) 


x4(t)  =  -(1000x2(t)+350x3(t)+35x4(t)+27  z^t)), 

and    the   output   of    the    system   is 

c    (t)    =  x..  (t)     .  (IV. 28) 

s  1 

2.  The  Feedback  Loop  (N2  -  Block) 

The  output  of  the  dead  zone  block  (see  Fig.  IV. 11(c)) 

can  be  represented  as 

0  if  |  xx(t)|  <   Kc 

z2(t)  =  ^x1(t)-Kc    if    xx(t)  >  Kc        (IV. 29) 

x1(t)+Kc    if    xx(t)  <-Kc  . 

3 .  Error 
The  error  is  given  by 

e(t)  =  r(t)  -  z2(t)  .  (IV. 30) 
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4.   Compensators 

The  cascade  compensators  have  transfer   functions 
of  the  form 

,   s  +  k .p. 
Gc.(s)  =F"   s  +  p.1   '    i=i'2  <IV-31> 

in  the  time  domain;   they  can  be  represented  by  the  confi- 
guration shown  in   Fig.  IV. 12.   For  the  first  compensator 

x6(t)  =  -p1  x6(t)  +  p1(l-l/k1)  e(t)  ,  (IV. 32) 

and  its  output  is 

Yl(t)  =  e(t)/k1  +  x6(t)  .  (IV. 33) 

Similarly,  for  the  second  compensator 

x5(t)  =  -  p2(x5(t)  -  (l-l/k2)  Yl(t))  (IV. 34) 

and  its  output  is 


y2(t)  =  Y]L(t)A2  +  x5(t)  .  (IV. 35) 


5.   N-l  -  Block 

The   output   of    the   saturating   amplifier   is    a   non- 
linear   function   of    the  magnitude   of   its    input    (see   Fig. 
IV.ll(b));       its    output  can  be    represented  as 

K  if      z..  (t)     >  K 

x  1  x 

zn  (t)    =    {   K  if      z, (t)     ^  K  (IV. 36) 

1  n  In 

K  y2(t)         otherwise 
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A  performance  index,  which  measures  the  difference  of 
the  outputs  of  the  nonlinear  and  ideal  systems,  for  a  given 
magnitude  of  a  step  input  (R) ,  is  defined  as 

ftf 

J(R,0?)  =  xi;L(t)  =  J  (xAt)    -  x?(t))  dt  ,       (IV. 37) 

o 

and  minimized  with  respect  to  the  free  system  parameters 
{01   =  (k,  k,,  p,  p,,)  )  and  the  optimum  parameter  values  (#*) 
are  found. 

Since  the  real  system  contains  nonlinearities  which 
are  amplitude  dependent,  the  performance  index  and  the  op- 
timum parameter  values  are  functions  of  the  magnitude  of 
the  step  input  (R) .   For  various  amplitudes  of  step  inputs 
the  optimum  parameter  values  are  tabulated  in  Table  IV. 2. 
The  time  responses  of  the  ideal  and  real  systems  for  the 
conditions  given  in  the  last  three  lines  of  Table  IV. 2  are 
shown  in  Figs.  IV. 13,  IV. 14  and  IV. 15. 

If  it  is  known  that  the  actual  system  is  going  to  be 
operated  usually  with  step  inputs  of  a  certain  amplitude 
(say,  R  =  2.5),  then  the  design  can  be  completed  by  using 
the  last  line  of  Table  IV. 2;   in  this  case  the  dynamic  re- 
sponse of  the  system  may  not  be  satisfactory  for  step  inputs 
much  larger  or  smaller  than  the  selected  value  (R  =  2.5). 
To  study  the  effect  of  the  input  magnitude,  the  pole  loca- 
tions are  fixed  at  the  optimum  value  found  for  the  step 
input  with  2.5  magnitude, 

p  =  10.2764  =  Constant  ,  (IV. 38) 
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STEP  INPUT 

MAGNITUDE (R) 

kl  '  k2 

Pl 

,         P2 

0.5 

1.7234 

8. 

.2500 

1.0 

1.2012 

8. 

7051 

2.0 

1.1340 

10. 

3790 

2.5 

1.1357 

10. 

2764 

TABLE    IV. 2 


STEP  INPUT 

MAGNITUDE 

Kl  '  K2 

0.1 

1.8115 

0.2 

1.8008 

0.3 

1.7754 

0.4 

1.7480 

0.5 

1.7215 

0.6 

1.6904 

0.7 

1.2139 

0.8 

1.2100 

0.9 

1.2021 

1.0 

1.1885 

1.5 

1.1502 

2.0 

1.1336 

2.5 

1.1355 

3.0 

1.1512 

3.5 

1.1678 

4.0 

1.3162 

4.5 

1.3543 

5.0 

1.3787 

TABLE    IV. 3 
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Fig.  IV. 13.  EXAMPLE  11.   Compensated  (A)  and  desired  (B)  sys^ 

terns'  responses  (r(t)  »  ll(t)). 
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Fig.  IV. 14.  EXAMPLE  11.   Compensated  (A)  and  the  desired  (B) 

systems'  responses  (r(t)  =  2  xll(t)) 
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and  the  optimum  compensation  ratios  (k*,k*)  are  found  and 
tabulated  (see  Table  IV. 3)  for  various  step  input  magnitu- 
des (R)  . 

A   graph  showing  the  relationship  between  the  step  in- 
put magnitudes  and  the  compensation  ratios  (R  vs.  (k,  and 
k2) )  is  plotted  (see  Fig.  IV. 16).   This  curve  can  be  divid- 
ed into  three  regions: 

(1)  0.1  <  R  <  0.6 

(2)  0.7  ^  R  ^  3.5  (IV. 39) 

(3)  3.5  <  R  ^  5. 

where  R  is  the  magnitude  of  the  step  input.   The  flat  sec- 
ond region  indicates  that  if  a  suitable  value  for  the  com- 
pensation ratios  is  chosen,  the  system  has  satisfactory 
output  responses  for  the  input  magnitudes  in  the  second 
region.   For  example, 

k,  =  k9  =  1.18  , 
1     Z  (IV. 40) 

Pl  =  P2  =  10*2764 

may  be  chosen  as  a  (fixed)  set  of  design  parameters  for  the 
input  magnitudes  in  Region  2.   If  this  range  covers  most  of 
the  operational  input  values  of  the  system  and  some  degrada- 
tion from  the  desired  response  can  be   tolerated   for  the 


1 

Since  the  system  has  a  stable  limit  cycle,  growing  oscil- 
lations do  not  occur  for  any  input  magnitude. 
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Input  magnitude  (R)  v.s.  Optimum 
compensation  ratios  (K*  and  K*) 
for  pole  locations  are  fixed 
at  P,  =  p2  =  10.264. 
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other  input  magnitudes,  this  design  is  acceptable;  other- 
wise another  compensation  scheme  (probably  adaptive  or  non- 
linear) must  be  considered. 

E.    CONCLUDING  REMARKS 

The  method  "Optimization  for  the  Best  Response  in  the 
Time-Domain"  is  directly  applicable  to  the  design  of  sampled 
data  systems  and  systems  with  one  or  more  nonlinear  elements 
or  time  delays. 
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V.   A  TIME-DOMAIN  STABILITY  ANALYSIS  METHOD 
USING  NUMERICAL  TECHNIQUES 

A.  GENERAL 

A  numerical  method  for  the  stability  analysis  of  feed- 
back control  systems  is  introduced  in  this  chapter.   This 
numerical  method  is  applicable  to  a  large  variety  of  feed- 
back control  systems,  but  it  may  be  valuable  especially 
for  the  nonlinear  and  sampled-data  systems  and  systems  with 
time  delay,  since  general  and  practical  methods  are  rare  in 
these  fields. 

If  the  stability  analysis  is  carried  out  for  a  system 
which  has  been  designed  by  using  the  method  "Optimization 
for  the  Best  Response  in  the  Time-Domain",  then  only  minor 
changes  in  the  computer  program  yield  the  stability  limit 
of  the  system;   for  this  reason,  this  method  is  especially 
convenient  for  these  types  of  design  problems. 

B.  GENERAL  PHILOSOPHY 

The  theory  behind  the  numerical  method  is  quite  simple 
It  again  depends  on  the  idea  of  "shaping  the  transient  re- 
sponse". The  state  equations  of  the  system  are  integrated 
with  respect  to  time  by  using  a  numerical  method,  the  out- 
put of  the  system  is  observed  during  the  integration  (from 
zero  time  to  final  time  -  t _ ) ,  and  the  maximum  and  minimum 
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points  of  the  output  are  found.  A  performance  index,  to 
measure  the  sum  of  the  differences  of  all  the  maxima  and 
the  differences  of  all  the  minima  is  defined 

n 

J(K)  =   ) (x  j   -  x  =>-1)2  +  (x  ?      -  x  J'1)2    ,     (V.l) 
v  '  Ls   max    max  '  v  mm    min  '-      '  v    ' 

j=2 

where  K  is  assumed  to  be  the  only  free  variable,  the  super- 
script j  is  the  counter  for  the  maxima  and  minima,  and  n 
is  the  number  of  maxima  and  minima  found  in  the  time  inter- 
val of  interest. 

Minimization  of  this  performance  index  with  respect  to 
the  free  variable  (K)  yields  the  value  of  K*,  which  puts 
the  system  at  its  stability  limit,  because  with  the  K*  value, 
the  output  of  the  system  becomes  constant-magnitude  oscil- 
lations. 

For  systems  with  two  variables,  a  stability  curve  can 
be  obtained  by  holding  one  of  the  variables  at  a  constant 
value  and  applying  the  above  procedure;   at  the  end  of  the 
minimization,  one  point  is  found  on  the  stability  curve  (in 
the  system's  parameter  plane).   By  assigning  other  values 
to  the  fixed  variable  and  minimizing  the  performance  index 
with  respect  to  the  free  variable,  many  points  of  the  sta- 
bility curve  can  be  obtained  and  plotted. 

The  same  idea  can  be  extended  to  systems  with  more 
than  two  variables  and  stability  surfaces  or  stability 
hypersurfaces  can  be  obtained,  but  this  may  require  too 
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much  computer  time,  and  the  practical  use  of  stability 
hypersurfaces  in  n-dimensional  Euclidean  spaces  is  ques- 
tionable.  For  practical  design  and  analysis  problems,  the 
intersection  of  the  stability  hypersurface  with  a  two-di- 
mensional plane  may  yield  the  desired  insight,  and  it  can 
easily  be  obtained. 

The  terms  in  the  performance  index  (values  of  the  suc- 
cessive maxima  and  minima)  are  found  by  adding  simple  logic 
to  the  subroutine  which  calculates  the  value  of  the  per- 
formance index.   The  maximum  and  minimum  points  of  the  re- 
sponse can  be  observed  by  different  methods  —  the  simplest 
logic  is  discussed  here.   At  the  end  of  each  integration 
step  (of  the  system's  state  equations)  the  slope  of  the 
system's  output  is  checked.   If  its  sign  is  altered  during 
the  integration  step,  or  if  it  takes  the  value  zero,  a  max- 
imum or  a  minimum  has  occurred.   For  a  continuous  output 
response  every  maximum  must  be  followed  by  a  minimum  and 
vice  versa;   for  this  reason  the  performance  index  defined 
in  Eq.  (IV. 1)  may  be  altered  to  the  simpler  form 

n 

J(K)  =   T(x-J  -  x-|-2)2  ,  (V.2) 

L->     m    m 

j=3 

where  x   represents  both  the  maxima  and  the  minima,  and  n 
m 

is  the  sum  of  the  number  of  maximum  and  minimum  points 
observed  in  the  time  interval  of  interest. 
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To  demonstrate  the  application  of  the  numerical  sta- 
bility analysis  method  three  examples  are  given. 

C.  EXAMPLE  12 

The  stability  limit  of  the  chemical  process  with  in- 
herent time  delay  is  to  be  found.   The  system  (Fig. IV. 1) 
and  its  state  equations  (IV. 2)  are  given  in  Example  9  and 
will  not  be  repeated  here. 

A  performance  index  as  defined  in  Eq.  (V.2)  is  mini- 
mized with  respect  to  the  only  free  system  parameter  (K 
the  amplifier  gain)  and  K*  =  2.125  is  found  for  the  sta- 
bility limit.   To  check  this  result,  the  system  is  simulat- 
ed with  this  gain  value  and  the  response  shown  in  Fig.  V.l 
is  obtained. 

D.  EXAMPLE  13 

A  sampled-data  system  (with  a  zero-order  hold)  was 
optimally  compensated  for  a  desired  time  response  by  using 
two  free  system  parameters  —  T  (the  sampling  period)  and 
K  (the  system  gain)  (see  Example  10).   The  system  is  shown 
in  Fig.  IV. 7 (a) . 

The  sampling  period  can  be  varied  from  0.1  second  to 
1.5  seconds  and  the  gain  can  be  varied  from  zero  to  15. 

The  value  of  the  sampling  period  is  varied  in  steps; 
it  is  first  set  to  0.1  then  increased  by  0.1  steps  (up  -o 
1.5) .   At  each  value  of  the  sampling  period,  a  performance 
index  similar  to  the  one  defined  in  Eq.  (V.2)  is  minimized 
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Fig.  V.l.  EXAMPLE  12.   Digital  simulation  of  the  system 

at  its  stability  limit  (K  =  2.125) 
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with  respect  to  K  and  the  stability  curve  in  the  T  -  K 
plane  is  plotted  (see  Fig.  V.2) .   This  curve  simply  gives 
the  maximum  allowable  gain  for  a  selected  value  of  the 
sampling  period,  or  vice  versa.   Several  points  on  the  sta- 
bility curve  are  checked  by  digital  simulation;  for  example, 
T  =  0.2  ,  K  =  10.34  are  read  from  Fig.  V.2  and  with  these 
parameter  values  the  system  is  simulated,  and  the  time  re- 
sponse shown  in  Fig.  V.3  is  obtained.   Similarly,  for 
T  =  1.5  and  K  =  1.74,  the  time  response  shown  in  Fig.  V.4 
is  obtained. 

E.    EXAMPLE  14 

The  stability  limit  of  the  system,  shown  in  Fig.V.5(a) 
is  to  be  studied.   The  uncompensated  linear  system  is  un- 
stable when  K  >  193.23  (see'  Example  2).   The  system  is  com- 
pensated by  using  a  tachometer,  and  K.  =  0.4  is  selected 
for  the  desired  response.   The  linear  compensated  system 
(with  K   =  0.4)  is  unstable  when  K  >  1020. 

The  tachometer  used  in  the  system  has  a  nonlinear 

characteristic  and  it  saturates  when  its  output  reaches  the 

saturation  level  (Z  ) .   In  other  words,  the  output  of  the 

m 

tachometer  is  a  function  of  the  magnitude  of  its  input  and 
can  be  represented  as  (see  Fig.  V.5(b)) 

(   Z  if    K.c  (t)  ^   Z 

m  t  s         m 

z(t)  =<-Z  if    K.c  (t)  <  -Z  (V.3) 

v  '         \       m  t  s        m 

K.c  (t)    otherwise 
V   t  sv 
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Fig.  V.3.  EXAMPLE  13.  Time  response  of  the  system  at 

one  of  its  stability  limits 
(T  =  0.2,  K  =  10.34)  . 
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Fig.  V.4.  EXAMPLE  13.  Time  response  of  the  system  at 

one  of  its  stability  limits 
(T  =  1.5  ,  K  =  1.74) . 
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For  a  qiven  saturation  level  (Z  ) ,  the  stability  of 
the  system  depends  on  the  value  of  the  gain  (K)  and  also  on 
the  magnitude  of  the  test  (input)  signal  (R) . 

The  state  equations  of  the  system,  written  directly 
from  Fig.  V.5(a),  are 


x1(t)  =  x2(t) 


x2(t)  =  x3(t) 
x3(t)  =  x4(t) 


(V.4) 


x4(t)  =  K(r(t)-  z(t)-  x1(t))-(50x2(t)+65x3(t)+16x4(t))/ 

where  the  input  is  a  step  function  with  magnitude  R,  i.e. 

r(t)  =  R  x  4(t)  .  (V.5) 

The  output  of  the  system  is 

cs(t)  =  xx(t)  ,  (V.6) 

and  its  time  derivative  is 

cs(t)  =  x2(t)  .  (V.7) 

For  a  given  step  input  magnitude  (R) ,  the  state  equa- 
tions of  the  system  are  integrated  by  using  a  numerical 
method;   after  each  integration  step,  the  logic  given  in 

(V.3)  is  used  to  calculate  z(t)  (Z   =  0.5  is  used).   The 

m 

sign  of  x~(t)  is  also  checked  against  its  previous  sign  and 
the  maximum  and  minimum  points  of  the  output  (x, (t))  are 
found.   The  performance  index  (Eq.  (V.2))  is  minimized  with 
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respect  to  (K) .   The  result  of  this  minimization  is  the 
maximum  allowable  gain  for  a  given  input  magnitude  (R)  and 
for  the  given  saturation  level  (Z   =  0.5).   After  a  minimi- 
zation is  completed,  another  magnitude  for  the  step  input  is 
chosen  and  the  minimization  is  repeated.   The  maximum  gain 
values  (K)  for  each  magnitude  of  the  step  input  (R)  are 
stored  and  the  stability  curve  in  the  R  -  K  plane  (see  Fig. 
V.6(a))  is  obtained  (for  Z   =0.5).   When  the  input  magnitu- 
de is  small  (R  ^  0.6),  the  operation  of  the  system  is  linear 
and  the  stability  limit  is  at  K  =  1020.   As  the  input  magni- 
tudes increase,  the  tachometer  feedback  becomes  less  effec- 
tive, due  to  saturation.   For  example,  when  the  step  input 
magnitude  (R)  is  100,  the  system  is  unstable  for  K  >  196.07; 
for  this  reason  the  stability  curve  (shown  in  Fig.  V.6(a)) 
starts  from  K  =  1020  for  small  input  magnitudes,  and  asymp- 
totically approaches  the  uncompensated  system's  stability 
limit  (K  =  193.23)  as  the  input  magnitudes  increase.   The 
same  stability  curve  of  Fig.  V.6(a)  is  replotted  to  show 
the  effect  of  magnitudes  up  to  R  =  100  (see  Fig.  V.6(b)). 

As  a  different  problem,  the  input  magnitude  is  fixed 
at 

R  =  2.5  =  Constant  ,  (V.8) 

and  the  effect  of  the  saturation  level  on  the  stability 

limit  is  studied. 

The  saturation  level  of  the  tachometer  feedback  (Z  ) 

m 

is  kept  constant,  and  the  minimization  is  carried  out  with 
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Fig.  V.6(a).  EXAMPLE  14.   The  stability  curve  in  the  "R  -  K' 

plane. 
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Fig.  V.6(b).  EXAMPLE  14.  The  stability  curve  in  the  extended 

"R  -  K"  plane. 
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respect  to  K;   when  the  value  of  K*  for  a  given  Z   is  found, 

3       m  .        ' 

Z   is  changed  to  another  value  and  the  minimization  is  re- 
in 

peated.   The  maximum  K  values  for  each  saturation  level  on 
the  tachometer  feedback,  in  the  range 

0.7  ^  Zm  ^  5.7  (V.9) 

m 

are  found  and  the  stability  curve  in  the  Z   -  K  plane  (for 
R  =  2.5)  is  plotted  in  Fig.  V.7. 

F.    CONCLUDING  REMARKS 

Two  words  of  caution  about  the  computer  programming 
are  appropriate  here.   The  performance  index  defined  in 
Eq.  (V.2)  has  at  least  two  minima;   one  of  these  occurs 
when  the  system  is  overdamped  and  the  time  response  has  no 
critical  points.   In  this  situation  both  the  value  of  the 
performance  index  and  its  gradient  become  zero,  and  this 
erroneously  indicates  a  point  on  the  stability  curve.  Some 
provision  should  be  added  to  the  computer  program  to  avoid 
this  erroneous  result.   Counting  the  number  of  maximum 
points  in  the  time  response,  and  finding  new  starting 
points  if  this  number  decreases  below  a  certain  predeter- 
mined constant,  is  one  way  of  avoiding  this  difficulty. 

Secondly,  the  performance  index  defined  in  Eq.  (V.2) 
measures  the  differences  in  successive  peaks  and  successive 
low  points  of  the  time  response;   but  it  is  also  a  function 
of  the  number  of  maxima  and  minima  which  occur  in  the  time 
interval  of  interest.   This  time  interval  may  contain 
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Fig.  V.7.  EXAMPLE  14.  The  stability  curve  in  the  "Z   -  K" 

plane.  m 
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different  numbers  of  maxima  and  minima  for  different  para- 
meter values.   To  make  the  performance  index  independent  of 
the  number  of  maxima  and  minima,  only  a  certain  number  of 
them  may  be  included  in  the  performance  index,  for  example 
the  third,  fourth  and  fifth  peaks  and  the  low  points  of  the 
time  response  may  be  measured  and  included  in  the  perform- 
ance index.   Averaging  the  performance  index  over  the  total 
number  of  maxima  and  minima  may  also  overcome  this  difficul- 
ty. 

The  two  programming  difficulties  discussed  above  are 
not  really  serious  and  with  a  little  care  they  can  be  over- 
come.  The  numerical  method  used  to  investigate  the  stabil- 
ity of  the  feedback  control  systems  can  be  very  attractive 
for  the  stability  analysis  of  nonlinear  and  sampled-data 
systems  and  systems  with  time  delay. 

This  method  is  an  extension  of  the  design  method 
"Optimization  for  the  Best  Response  in  the  Time-Domain". 
If  a  system  is  designed  by  using  this  method,  the  stabil- 
ity analysis  can  be  carried  out  with  little  additional  ef- 
fort, since  only  minor  changes  (primarily  in  the  perform- 
ance index)  convert  the  subroutine  which  was  used  in  op- 
timal design  to  the  one  to  be  used  in  the  stability  analy- 


sis 
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VI.   CONCLUSIONS 

Two  synthesis  methods  for  the  design  of  multi-parameter 
dynamic  systems  have  been  presented.   The  powerful  mathemat- 
ical tools  of  Optimal  Control  Theory  have  been  applied  to  a 
variety  of  system  design  problems. 

The  first  method,  "Optimization  for  the  Best  Root  Loca- 
tions  in  the  s-Domain"  has  been  applied  to  linear  system 
design  problems.   The  second  method,  "Optimization  for  the 
Best  Response  in  the  Time-Domain"  has  been  applied  to  the 
design  of  linear,  nonlinear  and  sampled  data  systems  and  to 
a  system  with  time  delay.   Eleven  design  examples  have  been 
considered;'   successful  solutions  of  these  examples  may  al- 
low one  to  conclude,  "These  design  methods  provide  fast, 
direct  and  efficient  solutions  to  practical  system  design 
problems,  and  they  can  compete  with  any  existing  design 
methods . " 

The  second  method  has  been  extended  to  a  numerical  sta- 
bility analysis  method,  and  its  effectiveness  has  been  de- 
monstrated by  three  examples. 

The  methods  presented  in  this  thesis  provide  computer- 
ized design  procedures  for  a  wide  class  of  dynamic  systems, 
and  hopefully,  they  provide  one  of  the  necessary  bridges  to 
narrow  the  existing  gap  between  the  theory  and  practice  of 
system  control. 
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13.  ABSTRAC  T 


Two  design  methods  for  multi-parameter  dynamic  systems  are  pro- 
posed.  They  are  intended  to  eliminate  the  limitations  and  disadvan- 
tages of  the  existing  design  methods.   The  powerful  mathematical  tools 
of  optimal  control  theory  are  applied  to  the  practical  design  problems 
of  classical  control. 

The  first  method  is  intended  for  linear  systems  only;   the  design 
problem  is  solved  in  the  s-domain,  by  finding  "the  best  root  loca- 
tions" of  the  system's  characteristic  equation.   In  the  second  method, 
the  design  problem  is  solved  by  finding  "the  best  response"  of  the  sys- 
tem in  the  time  domain.   The  second  method  is  applicable  to  a  wide 
range  of  dynamic  systems;   it  can  be  used  to  synthesize  linear,  non- 
linear and  sampled-data  systems,  and  systems  with  time  delay.   This 
method  is  also  extended  to  a  numerical  stability  analysis  procedure. 

Fourteen  examples  are  presented  to  illustrate  the  applications  of 
the  methods. 
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